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ABSTRACT 

Identified as extinction features against the brigfit Galactic mid-infrared back- 
ground, infrared dark clouds (IRDCs) are thought to harbor the very earliest 
stages of star and cluster formation. In order to better characterize the proper- 
ties of their embedded cores, we have obtained new 24 //m, 60-100 //m, and sub- 
millimeter continuum data toward a sample of 38 IRDCs. The 24 //m Spitzer 
images reveal that while the IRDCs remain dark, many of the cores are asso- 
ciated with bright 24 emission sources, which suggests that they contain 
one or more embedded protostars. Combining the 24 /xm, 60-100 /im, and sub- 
millimeter continuum data, we have constructed broadband spectral energy dis- 
tributions (SEDs) for 157 of the cores within these IRDCs and, using simple 
gray-body fits to the SEDs, have estimated their dust temperatures, emissiv- 
ities, opacities, bolometric luminosities, masses and densities. Based on their 
Spitzer /IKKC 3-8 /xm colors and the presence of 24 //m point source emission, 
we have separated cores that harbor active, high-mass star formation from cores 
that are quiescent. The active 'protostellar' cores typically have warmer dust 
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temperatures and higher bolo metric luminosities than the more quiescent, per- 
haps 'pre-protostellar', cores. Because the mass distributions of the populations 
are similar, however, we speculate that the active and quiescent cores may rep- 
resent different evolutionary stages of the same underlying population of cores. 
Although we cannot rule out low-mass star-formation in the quiescent cores, the 
most massive of them are excellent candidates for the 'high-mass starless core' 
phase, the very earliest in the formation of a high-mass star. 

Subject headings: dust, extinction - stars: formation - infrared: stars - submil- 
limeter 



1. Introduction 

The earliest phase of isolated low-mass star-formation occurs within Bok Globules. 
Viewed against background stars, Bok globules are identified as isolated, well-defined patches 
of optical obscuration and have typical visual extinctions. Ay, of 1-25 mag (Bok & Reilly 



1947). The individual precursor to a low-mass star, referred to as the 'pre-protostellar core', 
are found within Bok Globules. These pre-protostellar cores are typically small (~0.05pc) 
and dense (10^-10^ cm~^), with low temperatures (~10K) and low masses (0.5-5 Mq; e.g.. 



Myers & Benson 1983 Ward-Thompson et al. 1994). 



Because many low-mass star-forming regions are nearby, their pre-protostellar cores and 
protostars have been studied extensively. Moreover, these studies have the added benefit 
of achieving sufficiently high spatial resolutions to distinguish the individual protostars, 
allowing one to study and characterize their various evolutionary stages, i.e., Class-0, I, II, 



and III (Lada & Wilking 1984 Adams et al. 1987 Andre & Montmerle 1994). These stages 



are characterized by IR spectral energy distributions (SEDs) corresponding to increasing 
black-body temperatures. Moreover, differences in the shapes of their SEDs trace the gradual 
removal of circumstellar material surrounding the central protostar as it evolves and emerges 
from its natal core. 

While these early evolutionary stages are well-characterized for low-mass star-formation, 
it has been difficult to classify high-mass protostars in a similar manner. The combination 
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of large distances to high-mass star-forming regions, their rarity and rapid evolution, and 
the fact that most high-mass stars form deeply embedded in dense molecular clumps and 
within clusters with many lower-mass stars nearby, makes their identification and separation 
difficult. However, in a recent study combining MSX, IRAS, and sub-millimeter data toward 



42 regions of high-mass star formation, Molinari et al. (2008) have investigated the evolution 



of SEDs for young, high-mass protostars and has attempted to classify them via their SEDs. 
They find that objects in apparently different evolutionary stages occupy different areas in 
a bolometric luminosity versus envelope mass diagram, in a similar manner to the low-mass 
regime and, thus, conclude that high-mass star-formation may be a scaled up analog to 
low-mass star-formation. If this is the case, then perhaps the early stages of high-mass 
star-formation can also be characterized via differences in the SEDs of their dense molecular 
cores. 

In order to better characterize the properties of high-mass pre-protostellar and proto- 
stellar cores, a large sample of cores in the very earliest stages of high-mass star-formation is 
required. High-mass stars and clusters form from cold, dense molecular clumps within giant 



molecular clouds (Blitz 1991 Blitz & Williams 1999; Lada & Lada 2003). Recent studies 



suggest that the cold precursors of warm cluster-forming molecular clumps can be identified 
as 'infrared dark clouds' (IRDCs; Simon et al.||2006b Rathborne et al.|[2006 ). Surveys of the 
Galactic Plane at mid-IR wavelengths made with the ISO and MSX satellites first identified 
IRDCs as dark extinction features seen in absorption against the bright mid-IR emission 



arising from the Galactic background (Perault et al. 1996; Carey et al. 1998; Hennebelle 



et al.||2001 ). IRDCs are ubiquitous across the Galaxy (Simon et al. 2006a) and are typically 



long and very filamentary. They correspond to the densest parts of much larger giant molec- 



ular clouds (Simon et al. 2006b) and are characterized by high densities (> 10^ cm ^), high 
column densities (~10^^-10^^ cm~^), and low temperatures (< 25 K; Egan et aL]|l998 Carey 



et al. 1998, 2000). 



Recently, Rathborne et al. (2006) conducted a survey of the 1.2 mm continuum emission 



toward 38 IRDCs using MAMBO-II on the IRAM 30 m telescope. These IRDCs were selected 



from the catalog of MSX IRDC candidates (Simon et al. 2006a) and all have known kinematic 



distances, determined via the morphological match of ^^CO(l-O) emission to the mid-IR 
extinction ( Simon et al.||2006b Jackson et al.|[2006 ). In these IRDCs Rathborne et al. (2006) 
found 190 cores, 140 of which are cold, compact cores. These cold, compact cores have typical 
sizes of < 0.5 pc and masses of ~ 120 Mq (Rathborne et al. 2006). Indeed, millimeter and 
sub-millimeter observations toward other IRDCs suggest such compact cores are ubiquitous 
Lis fc Carlstrom||1994[ [Carey et al.]|2000| [Garay et al]|2004[ [Ormel et al" 



within IRDCs (e.g.. 



2005] [Beuther et al.||2005D . Because IRDCs are cold, their thermal dust emission will peak in 



the far-IR/sub-millimeter regime. At these wavelengths the dust emission is optically thin. 
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making this regime the best for probing their internal structure and reveahng their star- 
forming cores. If IRDCs are the high-mass analogue to Bok globules and the cold precursors 
to cluster-forming molecular clumps, then these dense cores may be the precursors to the 



stars (Rathborne et al. 2006). 



IRAM and JCMT molecular line spectra, Spitzer 3-8 /im and 24 fim continuum images, 
and GBT water and methanol maser spectra toward a sample of cores within IRDCs reveal 
that most contain little evidence for active star- formation, such cores are called 'quiescent' 



(Chambers et al. 2009). However, some of the IRDC cores do appear to be actively form- 
ing stars as they show broad molecular line emission, shocked gas, bright 24 /xm emission, 
and strong water and methanol maser emission (e.g., Rathborne et al. 2005; Wang et al. 



2006 Chambers et al. 2009). Each of these tracers provides independent evidence for star- 
formation, either indirectly (from the interaction between the protostar and the surrounding 
core traced through the broad line emission and the shocked gas) or more directly (from the 
heating of the dust surrounding the central protostar traced via the bright 24 /im emission). 
It is likely, therefore, that these particular cores contain protostars. Indeed, a number of 
protostars have already been identified within IRDCs and span a range in mass, from low- 



and intermediate- mass (Carey et al. 2000 Redman et al. 2003) to high-mass (Beuther et al. 



Rathborne et al. 


2005 


Pillai et al. 


2006 


Wang et al. 


2006 



In order to characterize the cores within IRDCs, we have conducted a large, multi- 
wavelength observational survey, combining IR, sub-millimeter, and millimeter continuum 
data. These data, when combined to make broadband SEDs, provide estimates of the core 
dust temperatures, dust emissivities, opacities, bolometric luminosities, and masses. Because 
the cores have SEDs that peak in the far-IR, to date, it has been very difficult to estimate 
bolometric luminosities of the individual protostars due to the uncertain extrapolation from 
much longer and/or shorter wavelengths. Indeed, many previous studies using MSX and 
IRAS fluxes have been forced to make assumptions about the relative contributions to the 



far-IR flux from individual objects within a star-forming region (e.g. Molinari et al. 2008). 
Thus, the inclusion of the sub-millimeter and far-IR data is critical to determine reliable 
SEDs for the cores and to better constrain these parameters. In this paper we provide the 
first mid-IR, far-IR, sub-millimeter, and millimeter SEDs for a large sample of cores within 
IRDCs. 



2. Observations 



The source list for these observations comprise 38 IRDCs and the 190 cores identified 
within them (Rathborne et al. 2006). To characterize the emission from these IRDCs and 



their cores, we have obtained continuum data at many wavelengths: 24 /im, 60-100 /xm, 
350 yum, 450 /im, 850 yum, and 1.2 mm|^ Table [l] gives a summary of the telescopes, instru- 
ments, wavelength range, observing dates, angular resolution and 1 a noise for the data. 



2.1. 24 /im continuum images 

The 24 /im continuum images were obtained using the Multi-band Imaging Photometer 



(MIPS; Rieke et al. 2004) array on-board the Spitzer Space Telescope. Images toward 30 of 
the IRDCs were obtained as part of a Spitzer cycle 1 General Observer (GO) proposal. The 
images were obtained in the raster-scanning mode during the periods 2004 October 16-19 
and 2005 April 7-13. Because the IRDCs have complex and extended morphologies, the 
map sizes varied to cover the extent of the extinction at 8 /im. For 25 of the IRDCs, a 3 
column X 3 row (13'xl3') map was sufficient to cover the 8 fim extinction. For four of the 
remaining IRDCs we used a 3 column x 5 row (13'xl9') map, while toward the remaining 
IRDC, we used a 5 column x 5 row (19'x 19') map to image the full extent of the 8 /xm 
extinction. In all cases, the raster map was stepped by half the array between consecutive 
scans. Three repeats of the map (~3 sec exposure time per point) were combined to produce 
the final image which achieved a 1 cr point source sensitivity of 124 /iJy and an extended 
source sensitivity of 0.13MJysr^^ (calculated using the on-line sensitivity calculator). 

For the 8 IRDCs that were not part of our GO Spitzer proposal, we use 24/im images 
from MIPSGAL (see Carey et al.||2005" for a description of the survey details). These images 



were obtained by mapping large portions of the Galactic Plane using the 'fast' scanning 
mode. The resulting mosaics have a 1 a point source sensitivity of 207 /iJy and an extended 
source sensitivity of 0.21MJysr~^. 

All data were processed using the S13.2.0 version of the MIPS reduction pipeline. 



2.2. 60—100 /im continuum spectra 

We used the SED mode of Spitzer/MlPS to obtain long-slit, low-resolution (R~15-25) 
spectra toward a sample of these IRDC cores. Due to the saturation limits of MIPS toward 
the Galactic plane, our observations were limited to include only those cores with 1.2 mm 
fiuxes < 2 Jy. In total, we obtained a far-IR continuum spectrum toward 72 cores. For 65 
of the cores, 4 repeats of a 10 sec integration were combined to produce the final spectrum 



•^For details of the 1.2 mm contimium emission data see 



Rathborne et al. 



(2006). 
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which achieved a la sensitivity of ~ 30, 70, and 150 mJy at 60, 75, and 90 /xm respectively. 
For the remaining 7 bright cores (those with 1.2 mm fluxes between 1-2 Jy), 4 repeats of a 
3 sec exposure were combined. All spectra were obtained in the pointed observation mode. 
At these wavelengths, Spitzer has an angular resolution of ~13-24" (9.8 "pixels). 

Because of the contamination by the second order diffracted light and an inoperative 
detector module, the wavelength coverage of the spectra was restricted to 65-97/im. The 
spectra were acquired in two epochs: 2006 October and 2007 May. Data from 2006 were 
processed using the S 14.4.0 pipeline version, while 2007 data were reduced using version 
S16.1.0. The pointed SED-modc observation provides a set of six pairs of data frames 
between the target position ('on") and nearby sky position ('off'). For all analysis, we use 
the pipeline produced post-basic calibrated data (post-BCDs), which deliver mosaic images 
for the 'on' and 'off' spectra. 

We set the scan mirror to chop between the 'on' and 'off' positions with a chop throw 
of V to 3'. The 'off' position was selected for each individual core to be nearby and free 
from 1.2 mm continuum emission. However, because the chop distance and visibility of the 
cores were limited, many of the 'off' positions fall within the 8 /im extinction associated 
with the IRDC. Although the cores dominate, we found that many of the 'off' positions 
were not completely free of far-IR emission. Thus, we could not use the pipeline-produced 
'on-off' spectra due to the 'off' contamination. Instead, we generated our own background 
and sky-subtracted spectrum using the local background immediately surrounding the core. 

The spectra were ffux calibrated using data obtained from multiple observations of three 
bright stars, with 70 /im flux densities of 13-19 Jy. For point sources, the fluxes are accurate 
to within ~ 20%. The typical 1 a noise in these data is ~ 100 mJy beam~^. 

2.3. 350 /xm, 450 /xm, and 850 /xm continuum images 

The sub-millimeter continuum images were obtained using both the Caltech Submil- 
hmeter Observatory (CSO) and the James Clerk Maxwell Telescope (JCMT) via either a 
'scan-mapping' or 'on-the-fly' mode. Due to weather and time constraints, not all 38 IRDCs 
were observed at all sub-millimeter wavelengths: in total, 37 of the IRDCs were observed 
at 350 /im, 29 at 450 /im and 8 at 850 /im. All 350 fim continuum images were obtained 
with the CSO and all 850 /im continuum images were obtained with the JCMT. The 450 /im 
continuum images were obtained at either the CSO or the JCMT. 

The CSO data were obtained over three observing runs, 2004 September 5 and 9, 2005 
April 12-15, and 2006 April 25-28 while the JCMT data were obtained on 2004 September 
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12-13. During the CSO observations, the r225GHz measured from the CSO sky dipper was 
~ 0.07 during the 2005 and 2006 April observing periods and shghtly higher, 0.09, for 2004 
September. The measured 1 cr noise in the images is ~ 0.2 Jy beam"^. The JCMT data 
were obtained in better conditions which resuhed in a 1 cr noise of ~0.06 Jy beam"^. 

Where possible, the size of the individual maps were selected to cover the majority of the 
extinction feature and the 1.2 mm continuum emission associated with the IRDC and ranged 
from 2'x2' to 8'x8'. Standard reduction methods within the software packages CRUSH (for 
the CSO data) and SURF (for the JCMT data) were used to correct for atmospheric opacity 
and to remove atmospheric fluctuations. The data were flux-calibrated using either G34.3 
or Uranus. Due to the uncertainties in accurate flux-calibration of sub-millimeter data, the 
calibration of the fluxes quoted here have errors of ~ 40%. 



3. Results 

3.1. Continuum images 

Figure [T] shows an example of the continuum data obtained toward one of the IRDCs: 
G028. 53— 00.25. Here we show the 24 /xm, 850 /im, 450 /im, and 350 /xm continuum images 
overlaid with contours of the 1.2 mm continuum emission. The continuum data toward this 
IRDC is typical in that it shows that IRDCs remain dark at 24 yum, which confirms their 
extremely high column densities and low temperatures (see the appendix for all the 24 /xm 
images). Moreover, the 24 /xm extinction morphologies closely match the morphologies of 
the 1.2 mm continuum emission. This is expected as the 1.2 mm continuum emission traces 
the cold, dense dust which is blocking the Galactic mid-IR emission. While the IRDCs 
themselves remain dark, some IRDCs appear to contain bright, compact 24 /xm emission 
sources. In many cases these 24 /im emission sources are coincident with the millimeter 



cores (e.g.. Figs. 



A-4, 


A-7, 


A-9, 


A-25 



Sub-millimeter continuum emission was detected toward all the observed IRDCs and 
matches well the morphology of the 1.2 mm continuum emission (e.g. Fig. |l|. While the 
extended sub- millimeter continuum emission from the larger IRDC is sometimes faint, the 
cores are associated with bright, compact sub-millimeter continuum emission. Because strong 
emission at millimeter/sub- millimeter wavelengths can trace either temperature or column 
density enhancements, millimeter /sub-millimeter data alone are unable to accurately indi- 
cate the physical properties of potential star-forming cores. Thus, one needs an additional 
method to distinguish high temperature cores from high column density cores. 



- 8 - 



3.2. Spectral Energy Distributions 

To characterize the emission arising from the cores, we have constructed broadband 
SEDs by combining the 24 /xm, 60-100 /xm, 350 /im, 450 /xm, 850 /im, and 1.2mm continuum 
data. To model their SEDs we assume that the emission is optically thin and use a single 
temperature, modified gray-body function of the form, 

= ^],S,(T^)(l-e-^-) 

where fig is the source solid angle, B^[T£,) is the Planck function at the dust temperature 



(Td), and T,y is the optical depth (Gordon 1995). We calculate T250 using the relation T,y 



r25o(250/A)'', where /3 is the dust emissivity index and A is the wavelength in /xm. The free 
parameters in the fits are the dust temperatures (To), optical depths at 250 /im (r25o), and 
dust emissivity indices, (3. 

While this model may not be appropriate to characterize the emission in all cases, it is 
the simplest and most systematic approach we can take to obtaining estimates of the global 
dust parameters. Because we are interested in the entire core which is typically larger than 
protostellar envelopes and disks, this model will be sufficient and will produce a good approx- 
imation to the dust parameters on ~ 1 pc spatial scales. Give the limited measurements at 
different wavelengths in the majority of the SEDs, complex radiation transfer modeling that 



includes disks, envelopes, protostars as well as several geometrical effects (e.g. Robitaille et 



al. 2007) would not produce robust results. Moreover, because the Robitaille et al. (2007) 
models currently do not incorporate the data at wavelengths > 100 fim accurately, the use 
of these detailed models is not appropriate for our purposes (priv. com. T. Robitaille). 

The millimeter/sub- millimeter fluxes and sizes for each core were determined by fltting 
a two-dimensional gaussian proflle with a constant background level to the emission at each 



core position (determined from the 1.2 mm continuum emission; see Table A-1 and Rathborne 
etaL][2006| ). Each core was flt individually and the gaussian was inspected to determine that 
it fit well the fiux from the core while excluding the underlying large-scale emission from the 
cloud. 

The 24 /im fiuxes were calculated by interactively fitting a gaussian profile with a con- 
stant background level to the radial profile of the point source emission at each core position. 
This fitting takes into account any local variations in the noise of the data or variations in 
the local background emission. We assume an error for the 24 /im fiux calibration of 20%. 
In the majority of cases where a gaussian profile could not be fit to the emission, we have de- 
termined an upper limit to the fiux by calculating the mean value in a small region centered 
at the core position. Because the emission may not be optically thin at 24 /im, the 24 fim 
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fluxes are only included in the SED fitting for cores that have no 60-100 /xm continuum spec- 
trum. Because we wish to constrain the total luminosity of the core including any embedded 
sources, and the emission from the cores typically peaks in the far-IR, the 60-100 fim data 
were crucial in determining a good fit, particularly in the cases where no 24 /zm emission was 
detected. For inclusion within the SED the 60-100 /im continuum spectra were separated 
into 4 bins. The bin centers are at 55, 65, 75, 85, and 95 //m and are 10 /xm wide. The 
mean fluxes within these ranges were used to produce the points that are included within 



the SEDs. Table A-1 lists the coordinates, distances, and fluxes measured for all 190 cores. 



A core SED was generated if fluxes could be extracted for at least two of the millimeter/ sub- 
millimeter wavelengths. Thus, not all cores have an SED and derived parameters. We use 
the least squares fltting routine in IDL, MPFITFUN, to determine the best flt to the data. 
This procedure has the advantage of allowing us to constrain the values for the input pa- 
rameters, include the errors in the fltting, and to correctly handle the upper limits at 24 /im. 
To model the gray-body emission we assume that the cores are isothermal and use a range 
in Td of 5-55 K, T250 of 0.0001-1.0, and (3 of 0.5-2.5. To account for the errors we assumed 
'standard' weighting for the data such that the weight, W = 1.0/error^, where the error 
includes both the calibration uncertainty and the image noise, added in quadrature. When 
the core had a 24 /im upper limit or a 60-100 /im continuum spectrum, the weighting for 
the 24 /im flux was set to 0.0 which forced the fltting routine to ignore this data point. In 
all of these cases the resulting flt was consistent with the upper limit (i.e. the function fell 
below the limit). 



Figures A-1 - A-32 show the 24 /im image (overlaid with the 1.2 mm continuum emission) 
for each IRDC as well as the broadband SEDs for the cores. In total, broadband SEDs were 
generated for 157 of the sample of 190 cores. The millimeter /sub-millimeter /far-IR fluxes 
are shown on these plots as flUed circles with the corresponding error bars. The 24 /im fluxes 
are plotted as either flUed circles (when included in the fltting), open circles (when excluded 
from the fltting), or as upper limit arrows. 

For the 100 cores that have millimeter/sub-millimeter emission and either 24 /im emis- 
sion or a 60-100 /im continuum spectrum, we plot the grey-body function (solid line) deter- 
mined from the best flt. Labeled on each plot are the IRDC and core name, the output To, 
r, (3, luminosity, mass and the returned from the flt. Using the output errors for the To, 
r, P, we also include the functions that correspond to the upper and lower values for the 



each of the three parameters (dotted lines). Table A-2 lists the parameters determined for 
these cores. 



The remaining 57 cores have an upper limit at 24 /im and no 60-100 /im continuum 
spectrum, i.e., only fluxes in the millimeter/sub- millimeter regime. Because the emission 
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from the cores peaks in the far-IR, the absence of either a 24 /im flux measurement or far- 
IR continuum spectrum made the determination of the individual core parameters difficult. 
For these cores the plots show two functions (examples are G015. 31-00. 16 MM4, MM5, 
G025.04-00.20 MM3, MM4, MM5; Figs. |A^ and [A^ respectively) . The first fit (solid 
line) is determined using only the millimeter /sub-millimeter fiuxes. Because these SEDs 
only contain two data points and the emission likely does not cover the peak, the parameters 
determined from them are probably lower limits. The second fit (dashed line) was determined 
by including the 24 fim upper limit as a real data point. In this case, the values determined 
from the fitting represent upper limits. For these cores we include two values for the Td, t, /3, 
luminosity, mass and the plots; the estimated lower and upper limits. Table A-3 lists 

the lower and upper limits to the parameters for these cores. Because these fits are unreliable 
the parameters derived from them are not included in any of the following analysis. 



Discussion 



4.1. Evidence for active star formation: separation of cores with high-mass 

protostars 

To identify cores that may contain active star formation, we use a combination of their 
Spitzer /IRAC 3-8 fim colors and the presence or absence of compact 24 /im emission. We 



use the IRAC data and classification scheme presented in Chambers et al. (2009) but use 



our own 24 /im fiux measurements which were determined via a gaussian fit to the radial 



profile of the emission (rather than via aperture photometry as in Chambers et al. 2009). 



The IRAC classification scheme is based on the colors of objects in the IRAC 3-8 fim 
images which are grouped in to three categories: 'red', 'green', and 'blue'. Most objects have 
either red or blue colors in these bands which indicate that that the fiux is either increasing 
or decreasing, respectively, with wavelength. The 'red' objects are associated with bright 
8 fim emission and likely correspond to Hii regions while the 'blue' objects are associated 
with a region of bright 3.6 fim emission and are predominantly unextincted stars. Those 
rare objects that are 'green' show enhanced emission at 4.5 fim and could correspond to 
either an extincted star or, if the emission is extended, shocked gas. These latter objects are 
referred to as either 'green fuzzies' or EGOs ( [Chambers et al.|[2009| [Cyganowski et"aL|[2008| ) , 
and are thought to trace the shocked gas in an outflow ( Marston et al.||2004 Noriega- Crespo 



et al. 2004). Of the 190 cores within our sample. Chambers et al. (2009) find that 35 are 



associated with 'red', 47 with 'green', 6 with 'blue' objects, while the remaining 102 cores 
are not associated with any significant IRAC emission. 
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Emission at 24 /iin is an indicator of active star formation. Because this emission traces 
warm dust that is the heated as material accretes from a core onto a central protostar, the 
detection of a bright 24 /im point source associated with a dense core suggests that star 
formation is occurring within it. We find that 93 of the 190 cores have a detectable 24 /xm 
point source emission. 

Thus, the combination of a green fuzzy and 24 /xm emission toward a dense core suggests 
that there is both shocked gas and an accreting protostar within the core. Conversely, cores 
that do not contain an accreting protostar will remain cold, with dust temperatures too low 
to emit detectable 24 /im emission. Because they also lack outflows, they will also show 
no evidence of shocked gas. As a result, the detection of a green fuzzy and a 24 /im point 
source coincident with a core indicates active star formation within that core. Conversely, 
the non-detection of shocked gas and 24 /im emission indicates that the core may be cold, 
with no (detectable) internal exciting source. While the absence of these star forming tracers 
is not definitive proof that there is no current star formation occurring within the core, it is 
suggestive. For example, low-mass protostars may lack sufficient luminosity to be detectable 
at 24 /im given the typical distances to these IRDCs. However, a lack of heated dust and no 
shocked gas is certainly what is expected for a core in the cold, pre-protostellar phase. 



Using their IRAC 3-8 /im colors and their 24 /im emission we use the Chambers et al. 



(2009) technique to separate the cores into five specific groups: quiescent cores, intermediate 
cores, active cores, red cores and blue cores. We define a 'quiescent core' to be a core that 
contains no significant 3-8 /xm nor 24 /im emission. 'Intermediate cores' contain either a 
green fuzzy or a 24 //.m point source, but not both. 'Active cores' contain a green fuzzy and 
a 24 /im point source. 'Red cores' are those cores associated with bright 8/im emission. 
Because bright emission at 8 /xm will be dominated by emission from either heated dust 
or UV-excited polycyclic aromatic hydrocarbons (PAHs), these cores likely correspond to 
embedded H ii regions where high-mass stars have already formed. Because the intermediate 
cores are associated with only one of the two proposed criteria for active star formation, we 
have separated these into a distinct group. These cores perhaps represent a transition phase 
between the active and quiescent cores. 

Of the 190 cores, we group 35 as 'red', 6 as 'blue', 38 as 'active', 32 as 'intermediate', 
79 as 'quiescent'. We generated SEDs for 25 of the 'red', 3 of the 'blue', 34 of 'active', 29 of 
the 'intermediate', and 66 of the 'quiescent' cores. The IRAC classifications, the presence of 
24 /im emission and the resulting core designations are included on the SED plots. 
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4.2. Dust temperatures, opacities, and dust emissivities 

Figures [2]-[4] show histograms for the number distributions of the derived T/), r25o, and 
/3 for our sample of cores. These figures are separated into four panels and show the number 
distributions for the red, active, intermediate, and quiescent cores for which reliable SEDs 
could be generatecQ For these 100 cores, 22 can be classified as 'red', 34 as 'active', 25 
as 'intermediate', and 17 as 'quiescent'. Table |2] lists the median and standard deviation 
measured from these distributions. 

The derivation of core dust temperatures allows us to establish whether significant 
heating, either by accretion or embedded sources, is occurring within the cores (mean radius 
of ~ 0.5 pc). As expected. Figure [2] shows that the highest dust temperatures are measured 
toward the red cores (median T/5 of 41 K), which presumably already contain a high-mass 
star. The measured dust temperatures are cooler for the samples of active, intermediate, 
and quiescent cores (median of 34 K, 31 K and 23 K respectively) as expected if the star 
formation activity is in an earlier state. 

While these dust temperatures are derived for each core individually, one must exercise 
caution when applying these temperatures to the cores where there is clearly active star 
formation. For instance, if the red cores really do harbor a high-mass star and Hii region, 
then close to the protostar their internal temperatures should be heated to well above 41 K. 
Because our angular resolution (~ 11"; 0.2 pc) samples dust from both the central proto- 
star and its surrounding larger core, the derived dust temperatures reflect a beam-weighted 
average of the actual dust temperature distribution. If a cold ~ 1 pc core contains a small 
volume of heated dust, the derived dust temperature will be larger than that of the cold 
envelope but smaller than that of the heated region. While the general trend of higher dust 
temperatures for cores with active star formation compared to lower temperatures for cores 
without star formation is expected, it is likely that, within the star-forming cores, there 
is a small volume of gas which is heated to temperatures significantly greater than those 
calculated here. Because the dust and gas are unlikely to be completely thermally coupled, 
it is possible that in these cases the gas temperatures are significantly higher than the dust 
temperatures derived from the SED fit. 

The derived values of t^^q reveal that all cores are optically thin at 250 /xm (r25o ~ 0.01) 
regardless of their star- formation activity (Fig. |3]). The derived values for /3 (Fig. |4]) show 
a range in values across the complete range input into the gray-body fitting routine. Since 



^Because there are only 2 blue cores and they are presumably associated with foreground stars we do not 
include them here. 
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active star formation may change the properties of the dust grains, one might expect the 
emissivity index to vary between the different groups. It appears that, ahhough they have 
fairly similar median values, the distributions may be slightly different for the red and active 
cores compared to the distribution for cores with no apparent star-formation. 



4.3. Core masses, column and volume densities 

Mass estimates obtained from molecular line emission toward very cold, very high col- 
umn density cores are often unreliable. The combination of high molecular line optical depths 
toward these cores and the potential for molecular depletion makes it difficult to accurately 
trace the internal structure of such cores using molecular line emission. Because the emission 
from cold, dense dust peaks at millimeter /sub- millimeter wavelengths and is optically thin, 
it is a superior tracer of the internal structure and the masses of these cold, dense regions. 
However, accurate mass estimates from dust continuum emission require a good measure- 
ment of the dust temperatures and emissivities. Most mass estimates from dust emission 
assume a single value for To and (3. Because a large sample of cores potentially spans a 
large range in evolutionary stages and environments, the true values of and P may vary 
considerably within the sample. To better determine the masses, therefore, one requires 
estimates of the dust temperatures and emissivities for each core. We have achieved this 
for our sample of cores and, thus, can now obtain a more accurate census of their masses, 
column and volume densities. 

To calculate masses we use 

M 



k^B^{Td) 

( Hildebrand||l983 ), where is the observed integrated source flux density, D is the distance. 



Ky is the dust opacity per gram of dust, and By{T£,) is the Planck function at the dust 
temperature. In all cases, we assume a gas-to-dust mass ratio of 100. We use the individual 
values for To and /3 derived from the SED for each core. To calculate the masses and densities 
we use the flux measured at 1.2 mm and Ki.2mm, of 1.0 cm^ g~^ ( Ossenkopf fc Henning|[T994 ) . 



The beam-averaged H2 column density, N(H2), was calculated using the expression 



By{TD)pimpKi,Vti, 

where is the observed peak source flux density, yU is the mean molecular weight (2.8), irtp 
is the mass of a proton, and VL[, is the beam solid angle. The volume-averaged H2 density, 
n(H2), was estimated using the dust masses and by assuming a volume of a sphere. The 



derived masses, column and volume densities for the cores are listed in Table I A- 2 



- 14 - 



Figures |5]-[7] show histograms of the number distributions for the masses, column and 
volume densities. A summary of the median and standard deviation values calculated from 
these distributions is given in Table [2j We find that the median masses are comparable be- 
tween the red, active, intermediate, and quiescent cores (Fig. |5]). To quantify the similarities 
between these distributions, we have also calculated, via a Kolmogorov-Smirnov (KS) test, 
the probability that these distributions arise from the same parent population. We find that 
the probability that the red and active cores are derived from the same parent is 6%, the 
probability that the red and quiescent cores are derived from the same parent is 84%, while 
the probability that the active and quiescent cores are derived from the same parent is 39%. 

At all wavelengths the 'active' cores that show evidence for star formation activity have 
significantly more luminous emission compared to those that are more quiescent. However, 
such a close correspondence between the mass distributions of these populations imphes that 
the bright emission observed toward the more active cores may simply arise from differences 
in their internal temperatures and not because they have higher densities. Indeed, we find no 
significant difference between the distribution for the derived column and volume densities 
for the cores that show star formation activity compared to those that are more quiescent 
(Figs. |6| and (T)). 



4.4. Core luminosities: identifying high-mass protostars and high-mass 

starless cores 

Bolometric luminosities provide information on a core's embedded young stellar objects 
and evolutionary state. The core luminosities were estimated by integrating the emission 
under the best fit gray-body curve and using the derived kinematic distance. The luminosities 



are included on the SED plot for each core and are listed in Table |A-2 



Figure |8] shows the number distributions of the derived bolometric luminosities for the 
cores (the median and standard deviation values calculated from these distributions are given 
in Table [2]). While the complete sample of cores span a range in bolometric luminosities of 
~ 10-10^ L0, there is a clear trend for the core luminosity to decrease from the red cores 
to the quiescent cores. We find that the red cores have typical luminosities of ~ 10^ '' Li^ 



0) 



while the active and intermediate cores have lower luminosities of ~ 10^'^ L(7) and ~ 10^'^ Lp 



'^'-'^^ '-^ -^0 
1.9 



respectively. Lower still, the quiescent cores have typical luminosities of ~ 10 Lq. 

A rough approximation of the final stellar mass can be made from the estimates of the 
core luminosity. Low-mass (M<2 Mq) stars never achieve luminosities > 100 Lq in their 



pre-main-sequence evolution (e.g. Palla & Stabler 1990). Because L > 100 Lq for the 
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majority of the cores, it is unlikely that, in these cases, a single low-mass star dominates 
the core luminosity. Given our sensitivity and angular resolution, we cannot rule out the 
possibility, however, that many of the cores may comprise a cluster of unresolved low-mass 
stars that, together, produce a larger core luminosity. 

On the other hand, high luminosities may arise because the core contains a high-mass 
protostar (M>8 Mq). Because the luminosity for a high-mass protostar (~10^-10^ Lq) is 
much larger than the luminosity from a low-mass protostar (~ 100 Lq), cores with lumi- 
nosities > 10^ Lq probably harbor a high-mass protostar. Since high-mass stars typically 



form in clusters surrounded by many lower mass stars (Lada & Lada 2003), lower mass 



protostars may also be forming within these cores and contributing to the overall observed 
luminosities. Regardless of the exact number of stars within the core, it is likely that the 
observed bolometric luminosity is dominated by the most massive protostar. 

Using a bolometric luminosity of 10"^ Lq as a rough threshold, we can identify cores that 
have sufficient luminosities to harbor a high-mass protostar. When applying this threshold 
to the sample of 'active' and 'intermediate' cores, we find that 6 cores have luminosities 
> 10^ Lq and, thus, may contain a high-mass protostar (five from the active core sample 
and one from the intermediate core sample). Indeed, high- angular resolution millimeter/sub- 
millimeter interferometry images toward two of these cores (G024. 33+00. 11 MMl and G034. 43+00. 24 MM; 



Figs. A-9 and A-25 respectively) show that both cores contain bright, compact structures 



that remain unresolved at ~ 0.03 pc angular resolution ( Rathborne et al.]|2007 2008). More- 
over, their high-angular resolution spectra show many complex molecular emission features 
indicative of hot molecular cores (HMCs); an early stage in the formation of individual 



high-mass stars (e.g. Garay & Lizano 1999). 



The majority of the remaining active and intermediate cores (45 cores; 76%) have lumi- 
nosities between 10^ and 10^ Lq. If these cores contain a single protostar or main sequence 
star, then they likely correspond to an intermediate mass star (2 < M < 8 Mq). The ac- 
tive/intermediate cores that have lower bolometric luminosities (i.e. L < 10^ Lq) may either 
be in an earlier evolutionary phase or be only forming low-mass protostars. They contain 
heated dust, as evidenced by the 24 /im emission, but the current luminosity appears too 
low for a high-mass protostar. In these cases, the observed luminosity may arise from the ac- 
cretion of material from the core onto the low-mass protostars. Indeed, recent work suggests 
low/intermediate mass protostars may continue to accrete material from their surroundings 



and grow to become high-mass protostars at the end of their evolution (e.g. Beuther & 



Sridharan 2007). Thus, we cannot rule out the possibility that these cores may eventually 



give rise to high-mass stars. 



As expected, the quiescent cores all have lower (L<10'^ Lq) bolometric luminosities. 
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presumably because there is no high- mass central source to heat the dust. While the number 
of quiescent cores included within these histograms is small (17), there are 62 additional cores 
within these IRDCs that also meet our criteria for a lack of active, high-mass star formation. 
These cores have no detectable 3-8 /zm nor 24 /xm emission but do not have the far-IR 
continuum data which is invaluable for determining their properties. As a result, they are 
not included in the analysis here, but are also likely pre-protostellar in nature. 

The quiescent cores have no apparent evidence for star-formation in the Spitzer data. 
One possibility for this absence of star-formation activity is that they are young but will 
eventually form stars. On the other hand, it may be that these cores are currently forming 
only low-mass stars that are not bright enough to detect. Thus, we cannot rule out the pos- 
sibility that the quiescent cores are currently forming low-mass protostars. Indeed, using the 



expression from the recent work of Krumholz & McKee (2008) that predicts the luminosity 



to mass ratio, L/M, of a core powered by the accretion on to low mass stars, 

L/M = 3.6M-°-33s°-^6Tfe°-^6 Lq/Mq 

where M is the core mass in units of 100 Mq, S is the mass column density in gcm~^, and 
Tb is the background temperature in units of 10 K, we find that for the quiescent cores, 
the median value predicted for the luminosity due to the accretion on to low-mass stars is 
~300 Lq. There are 6 quiescent cores within our sample that have luminosities comparable 
to, or greater than, 300 Lq. Thus, their observed luminosity could be powered by accretion 
onto low-mass stars. 

At the distances to these IRDCs (typical distance of ~ 4 kpc) low-mass protostars are 
difficult to detect. Thus, many of the 'quiescent' cores may contain low-mass protostars that 
will not be detectable with the available sensitivity and angular resolution. If, on the other 
hand, these candidate quiescent cores are in fact in a pre-protostellar phase, the most massive 
cores are excellent examples for the elusive 'high- mass starless cores'; the very earliest stage 
in the formation of a high-mass star. 

For the analysis presented in this paper, we assume that each core contains a single 
protostar. Without sensitive, high-angular resolution data, it is difficult to determine any 
potential multiplicity within the cores. High-angular resolution (~ 2"; ~ 0.03 pc) millime- 
ter and sub-millimeter data toward 6 of the highest mass cores within the current sample 



reveal that 4 of the cores are resolved into multiple protostellar condensations (Rathborne 



et al. 2007, 2008). It may be that these cores will in fact give rise to star clusters. In all 
cases, however, it appears that one central, massive condensation most likely dominates the 
bolometric luminosity. 

Using the expression that relates the maximum stellar mass in a cluster {rrimax) to the 
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mass of the cluster (Mduster), ^max = 1-2 ^cz'-ufter (Larson 2003), we can estimate the mass 
of a cluster than will give rise to at least one high-mass star. Setting rrimax to 8 Mq, we 
find that the total cluster mass that will harbor a high-mass star should have a mass of 



~ 67 Mq. Assuming a cluster star formation efficiency of ~ 30% (Lada & Lada 2003), we 
calculate that such a stellar cluster could form from a core with a mass of ~ 225 Mq. Thus, 
in the pre-protostellar phase, a core that will form a star cluster that will eventually form 
a high-mass star ought to have a mass of at least ~ 225 Mq. Five quiescent cores within 
these IRDCs meet this criteria and, thus, will probably form a cluster with a high-mass star. 



In particular, the most massive of these, the cores G028. 53— 00.25 MMl and MM3 (Fig. A- 
7| have low dust temperatures (~ 20 K) and bolometric luminosities (~ 10^ Lq), but high 



masses (~ 10'^ ^o); column densities (~ 6x10^'^ cm ^) and volume densities (~ 10^ cm ^). 
These values are consistent with the expected properties of high-mass starless cores. 

Recent theoretical work suggests that to form a high-mass star from a core the cloud's 
critical mass column density, E, should be > lgcm~^ and its luminosity-to-mass ratio, L/M, 
> 10 Lq/ Mq (Krumholz & McKee 2008). Under these conditions, further fragmentation 



is suppressed by accretion onto nearby lower-mass protostars and the subsequent heating 
of the surrounding material. It is thought that an increase in temperature to ~ 100 K is 
sufficient to increase the Jeans mass and, thus, halt the fragmentation process. Figure |9] 
plots L/M versus S for our sample of cores, with these thresholds overlaid. 

We find that 6 red cores and 3 active cores have L/M > 10 Lq/ Mq and S > lgcm~^ 
and, thus, have sufficient luminosity and mass to form a high-mass protostar (according to 



the Krumholz fc McKee||2008 criteria). For the remaining cores, it appears that the fragmen- 
tation and the formation of lower-mass protostars is still possible. If these conditions hold, 
then these cores may be resolved into multiple protostellar condensations in high-angular 
resolution data. Because the size scales of these cores are ~ 0.5 pc and are comparable to the 
size scales over which a star cluster may form rather than the individual stars, it is unclear 
if these criteria are applicable. However, their low S is consistent with the fact that they 
will likely fragment further. 



4.5. Evolutionary sequence 

The combination of data sets spanning wavelengths from the millimeter through the 
sub-millimeter and far-IR allows us to investigate the differences between the cores that 
show evidence for high-mass star formation activity from those that are more quiescent. 
Because we can derive both masses and temperatures from the data, we can separate their 
effects. One hypothetical possibility is that the core properties merely reflect an evolutionary 
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sequence. If this idea is correct, then the colder 'quiescent' cores are pre-protostellar and will 
evolve into the warmer, 'active' protostellar cores as protostars form within them and begin 
to heat them internally. In this case, only the temperature of the cores should vary, since 
the mass for the entire cluster-forming core is unlikely to change very much as it evolves. 
Thus, one would expect the active cores to be warmer than the quiescent cores, but that 
their mass distributions should be the same. 

An alternative possibility is that the cores that show evidence for star-formation simply 
contain more luminous, high-mass protostars, whereas the more quiescent cores contain only 
lower mass protostars. Thus, the lack of obvious star-formation activity in the quiescent 
cores would then result from the difficulty to detect the much fainter lower-mass protostars. 
Since low-mass stars form from lower-mass cores, in this scenario, one would expect the 
quiescent cores to have both lower temperatures and smaller masses on average than the 
active cores. Thus, the data allow us to distinguish between these two competing ideas: 
evolution and mass. 

While the derived dust temperatures and luminosities are higher for the cores that show 
evidence for star-formation compared to those that are more quiescent, the mass distribu- 
tions for the two populations are similar. In fact, the KS test shows that the probability that 
the active and quiescent cores are derived from the same parent population is ~40%. The 
similarity in their mass distributions suggests that the underlying populations are nearly 
the same, and that the differences in their temperature and luminosity reflect different evo- 
lutionary stages. If this is the case, then the quiescent cores have enough mass to form a 
high-mass protostar, however, they simply have not yet begun the process. Such cores should 
still be cold, and since no high-mass protostars have formed yet, also have low bolometric 
luminosities. In contrast, the active cores contain an internal protostellar heating source(s) 
which results in their higher dust temperatures and luminosities. 

The evolutionary differences between the samples of cores can also be seen in a bolomet- 



10 


Sridharan et al. 


2002 


Molinari 



et al. 2008). For both the low- and high- mass regimes, sources in different evolutionary 
phases have been shown to lie in distinct regions within this diagram. As the protostellar 
activity increases, the luminosity also increases such that the different stages within the star 
formation process occupy regions that overlap in mass but are offset to higher luminosities. 



Similar to the results of Molinari et al. (2008), we find that these IRDC cores lie in the 



high-mass regime and that, for a given mass, the bolometric luminosities increase from the 
quiescent, to the active cores, to the red cores. This is consistent with the idea that the 
quiescent cores are in an earlier evolutionary phase compared to the active and red cores. 
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5. Conclusions 

To characterize the physical properties of cores within IRDCs we have obtained new 
24 /xm, 60-100 fira, and sub-millimeter continuum data toward a sample of 38 IRDCs. These 
IRDCs contain 190 compact cores, 140 of which are dark at 8 /xm and are cold, compact and 
dense. The Spitzer/MlPS 24 //m images reveal that while the IRDCs remain dark, many of 
their cores are associated with bright 24 /im emission sources. The sub-millimeter contin- 
uum data elucidate both the large- and small-scale structure of the IRDCs. Because emission 
at millimeter/sub- millimeter wavelengths can trace either temperature and/or density en- 
hancements one needs an additional method to distinguish between these two parameters. 
Using the presence or absence of 24 /xm point source emission, in combination with their 
Spitzer/IKAC 3-8 fim colors, we have classified the cores into five groups: red, active, inter- 
mediate, quiescent, and blue. We find that, of the 190 cores, 35 can be classified as red, 38 
as active, 32 as intermediate, 79 as quiescent, and 6 as blue. 

From gray-body fits to their spectral energy distributions (SEDs) we have determined 
the dust temperatures, emissivities, opacities, bolometric luminosities, and masses for a large 
sample of the IRDC cores. The derived distributions of the dust temperatures, luminosities, 
and masses for the different groups of cores reveals that the dust temperatures and lumi- 
nosities are higher for those cores that show active, high-mass star formation compared to 
those cores that are more quiescent. Lower dust temperatures and luminosities are expected 
for the quiescent cores because they presumably have no high-mass internal source to sig- 
nificantly heat the dust. Comparing the derived masses for the core samples, however, we 
find that their mass distributions are similar. We interpret this similarity to be a result of 
evolutionary differences: the cooler quiescent cores may be the pre-protostellar precursors 
to the warmer, more active protostellar cores. 

Using their derived bolometric luminosities, we estimate that ~10% of the cores that 
show evidence for star formation may contain high-mass protostars. If the quiescent cores are 
indeed devoid of star formation, then the most massive of these are excellent candidates for 
the 'high-mass starless core' phase, a very early phase in the formation of a high-mass star. 
Because of their distances, we cannot yet rule out the possibility that many of these cores 
may contain low-mass stars. Observations with ALMA will be crucial to address this issue. 
Nevertheless, our study supports the idea that IRDCs harbor the very earliest evolutionary 
stages in the formation of high-mass stars and, thus, clusters. 
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Appendix 



Figures A-1 A-32 show the 24 /zm image toward the IRDCs overlaid with the 1.2 mm 
continuum emission. Also included in these figures are the SEDs and gray-body fits for each 
of the cores that had sufficient data. Table A-1 list the coordinates, distances, and fluxes 
measured toward all cores. Marked on the SEDs and listed in Table IA-21 and Table IA-31 are 
the derived parameters. 
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Table 1: Summary of the observations. 



Telescope 


Instrument 


Wavelength 


Date 


Angular 


1 a 










Resolution 


sensitivity 


Spitzer 0.85 m 


MIPS 


24 /im 


2004 Oct, 2005 Apr 


6" 


~150 AtJy 


Spitzer 0.85 m 


SED-mode 


60-100 /^m 


2006 Oct, 2007 May 


13-24" 


~100 mjy 


CSO 10 m 


SHARC-2 


350, 450 Atm 


2005 Sept, 2005 Apr, 2006 Apr 


8" 


~ 200 mJy 


JCMT 15 m 


SCUBA 


450, 850 Atm 


2004 Sept 


8, 15" 


~ 60 mJy 


IRAM 30 m'' 


MAMBO-II 


1.2 mm 


2004 Feb 


11" 


~ 10 mJy 



"These data were presented in Rathborne et al. (2006). 



Table 2: Summary of the derived core properties 



Property Red Active Intermediate Quiescent 

Median Stddcv Median Stddev Median Stddev Median Stddev 



To (K) 


40.4 


5.7 


34.5 


6.5 


30.4 


8.1 


23.7 


5.3 


Log[T250] 


-1.894 


0.575 


-2.081 


0.550 


-2.456 


0.628 


-2.404 


0.501 


/3 


1.7 


0.4 


1.6 


0.4 


1.4 


0.5 


1.2 


0.5 


Log [Mass] (M0) 


1.85 


0.59 


2.06 


0.54 


1.97 


0.61 


1.92 


0.55 


Log[Luni] (Lg) 


3.65 


0.72 


2.88 


0.69 


2.27 


0.79 


1.84 


0.44 


E (gcrn^'^) 


-0.39 


0.39 


-0.15 


0.34 


-0.29 


0.42 


-0.25 


0.40 


Log[N(H2)] (cm-^) 


22.05 


0.42 


22.10 


0.43 


22.03 


0.36 


22.01 


0.29 


Log[n(H2)] (cm-a) 


5.82 


0.38 


5.98 


0.32 


5.81 


0.42 


6.06 


0.39 
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Fig. 1. — An example of the continuum data obtained toward the IRDC G028. 53— 00.25. 
From top left to bottom right the continuum images are: 24 /im, 350 /xm, 450 /xm, and 850 /xm. 
In all cases the contours overlaid are the 1.2 mm continuum emission from IRathborne et al.l 
( |2006[ ). The contour levels are 30 (3(t), 60, 90, 120, 180, 240 mJy beam^^. These images 
reveal that while the IRDC remains dark at 24 fim it is associated with bright emission at 
all sub-millimeter wavelengths. Moreover, the cores are clearly identified above the emission 
from the larger cloud. Broadband SEDs were generated for all the cores and their proper- 
ties determined from gray-body fits. See the appendix for the SED plots and the derived 
parameters. 



- 26 - 




Fig. 2. — Histograms of the dust temperature (Td) derived from the gray-body fits to the 
core SEDs. The panels shows the histograms for the samples of the red, active, intermediate, 
and quiescent cores. We find that, as expected, the derived dust temperatures decrease from 
the red, to the active, to the intermediate, to the quiescent cores. Table |2] lists the median 
and standard deviations of these distributions. 



-3.5 -3.0 -2.5 -2.0 -1.5 -1.0 -0.5 

Log[T25o] 



Fig. 3. — Histograms of the 250 fim opacity (r25o) derived from the gray-body fits to the core 
SEDs. The panels show the histograms for the samples of the red, active, intermediate, and 
quiescent cores. We find that the majority of the cores are optically thin, with low values 
of r25o (most have T250 < 0.01). Table |2] lists the median and standard deviations of these 
distributions. 
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Fig. 4. — Histograms of the dust emissivity index derived from the gray-body fits to the 
core SEDs. The panels show the histograms for the samples of the red, active, intermediate, 
and quiescent cores. Although they span a similar range, the median of the distribution of 
P for star-forming cores may be different from the median of the distribution of /3 for cores 
with no apparent star-formation . Table [2] lists the median and standard deviations of these 
distributions. 
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Fig. 5. — Histograms of the derived dust mass. The panels show the histograms for the 
samples of the red, active, intermediate, and quiescent cores. We find that the mass dis- 
tributions for the cores span a similar range, regardless of their apparent star-formation 
activity. Table [2] lists the median and standard deviations of these distributions. 
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Fig. 6. — Histograms of the derived H2 column density, N(H2). The panels show the his- 
tograms for the samples of the red, active, intermediate, and quiescent cores. We find that 
the cores have similar column densities, regardless of their apparent star-formation activity. 
Table [2] lists the median and standard deviations of these distributions. 
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Fig. 7. — Histograms of the derived H2 volume density, n(H2). The panels show the his- 
tograms for the samples of the red, active, intermediate, and quiescent cores. We find that 
the cores have similar volume densities, regardless of their apparent star-formation activity. 
Table [2] lists the median and standard deviations of these distributions. 
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Fig. 8. — Histograms of the bolometric luminosity derived from the gray-body fits to the 
core SEDs. The panels show the histograms for the samples of the red, active, intermediate, 
and quiescent cores. We find that, as expected, the derived luminosities decrease from the 
red, to the active, to the intermediate, to the quiescent cores. Table [2] lists the median and 
standard deviations of these distributions. 
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Fig. 9. — Plot of bolometric luminosity to mass ratio (L/M) versus mass column density 
(S) for the cores. The dotted lines mark the critical values of these parameters for the 



formation of a high-mass star (determined from the recent theoretical work of Krumholz & 



McKeepOOS ). These criteria mark the thresholds above which the fragmentation of a core is 



suppressed, by the increase in temperature and Jeans mass. According to these criteria, we 
find that nine of the cores may give rise to a high-mass protostar. For the remaining cores, it 
appears that the fragmentation and the formation of lower-mass protostars is still possible. 
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Fig. 10. — Plot of bolometric luminosity versus mass for the cores. We find that for a 
given mass the luminosity increases from the quiescent, to the active, to the red cores. This 
is consistent with the idea that the quiescent cores are in an earlier evolutionary phase 
compared to the active and red cores. 
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Fig. A-1. — Spitzer 24 yum image overlaid with 1.2 mm continuum emission for 
G015. 05+00. 07 (contour levels are 30 and 60mJy beam^^). The lower panels show the 
broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub- 
millimeter, and far-IR continuum data are shown as filled circles (with the corresponding 
error bars), while the 24 /im fluxes are shown as either a filled circle (when included within 
the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For cores that 
have measured fluxes only in the millimeter/sub- millimeter regime (i.e. a limit at 24 /im), 
we show the results from two fits: one using only the measured fluxes (solid line; lower limit), 
while the other includes the 24 yum limit as a real data (dashed line; upper limit). In all 
other cases, the solid line is the best fit gray-body, while the dotted lines correspond to the 
functions determined using the errors for the T^, r, and /3 output from the fitting. Labeled 
on each plot is the IRDC and core name, classification, and the derived parameters. 
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Fig. A-2. — Spitzer 24 /im image overlaid with 1.2 mm continuum emission for 
G015. 31— 00.16 (contour level is 30mJy beam^^). The lower panels show the broadband 
SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub-millimeter, 
and far-IR continuum data are shown as filled circles (with the corresponding error bars), 
while the 24 //m fluxes are shown as either a flUed circle (when included within the flt), an 
open circle (when excluded from the flt), or as an upper limit arrow. For cores that have 
measured fluxes only in the millimeter/sub-millimeter regime (i.e. a limit at 24 /xm), we 
show the results from two flts: one using only the measured fluxes (solid line; lower limit), 
while the other includes the 24 yum limit as a real data (dashed line; upper limit). In all 
other cases, the solid line is the best fit gray-body, while the dotted lines correspond to the 
functions determined using the errors for the T^, r, and /3 output from the fitting. Labeled 
on each plot is the IRDC and core name, classification, and the derived parameters. 
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Fig. A-3. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G018.82-00.28 (contour levels are 30, 60, 90, 120, 240 mJy beam-^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-IR continuum data arc shown as filled circles (with the 
corresponding error bars), while the 24 /im fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two flts: one using only the measured fluxes 
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Fig. A-4. — Spitzer 24 /im image overlaid with 1.2 mm continuum emission for 
G019.27+00.07 (contour levels are 30, 60, 90, 120, 240 mJy beam'^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /im fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter regime 
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Fig. A-5. — Spitzer 24 nm image overlaid with 1.2 mm continuum emission for 
G022.35+00.41 (contour levels are 30, 60, 90, 120, 240 mJy beam^i). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-IR continuum data are shown as fllled circles (with the 
corresponding error bars), while the 24 /xm fluxes are shown as either a fllled circle (when 
included within the flt), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fiuxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /im), we show the results from two fits: one using only the measured fiuxes 
(solid line; lower limit), while the other includes the 24 /im limit as a real data (dashed line; 
upper limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the T^, r, and /3 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A-6. — Spitzer 24 /im image overlaid with 1.2 mm continuum emission for 
G022. 73+00. 11 (contour levels are 20, 30 and 45mJy beam"^). The lower panels show 
the broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, 
sub-millimeter, and far-IR continuum data are shown as filled circles (with the correspond- 
ing error bars), while the 24 /im fluxes are shown as either a filled circle (when included 
within the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For 
cores that have measured fluxes only in the millimeter /sub- millimeter regime (i.e. a limit 
at 24 /im), we show the results from two fits: one using only the measured fluxes (solid line; 
lower limit), while the other includes the 24 /im limit as a real data (dashed line; upper 
limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
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Fig. A-7. — Spitzer 24 /.tm image overlaid with 1.2 mm continuum emission for 
G023.60+00.00 (contour levels are 30, 60, 90, 120, 240, 360, 480 mJy beam-^). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
milhmeter, sub-milhmeter, and far-lR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter /sub-millimeter regime 
(i.e. a hmit at 24 /xm), we show the results from two fits: one using only the measured fluxes 
(solid line; lower limit), while the other includes the 24 //m limit as a real data (dashed line; 
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Fig. A-8. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G024.08+00.04 (contour levels are 30, 60, 90, 120, 240 mJy beam"^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimctcr, and far-lR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /xm fluxes are shown as cither a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fiuxes 
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Fig. A-9. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G024.33+00.11 (contour levels arc 30, 60, 90, 120, 240, 360, 480, 840, 1200 mJy beam"^). 
The lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived 
from the millimeter, sub-millimeter, and far-IR continuum data arc shown as filled circles 
(with the corresponding error bars), while the 24 fim fluxes are shown as either a filled circle 
(when included within the fit), an open circle (when excluded from the fit), or as an up- 
per limit arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter 
regime (i.e. a limit at 24 /xm), we show the results from two fits: one using only the mea- 
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Fig. A-10. — Spitzer 24 //m image overlaid with 1.2 mm continmim emission for 
G024.60+00.08 (contour levels are 30, 60, 90, 120, 240, 360, 480, 840, 1200 mJy beam"^). 
The lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived 
from the millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles 
(with the corresponding error bars) , while the 24 /xm fluxes are shown as either a flUed circle 
(when included within the flt), an open circle (when excluded from the flt), or as an up- 
per limit arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter 
regime (i.e. a limit at 24 //m), we show the results from two flts: one using only the mea- 
sured fluxes (solid line; lower limit), while the other includes the 24 /xm limit as a real data 
(dashed line; upper limit). In all other cases, the solid line is the best flt gray-body, while 
the dotted lines correspond to the functions determined using the errors for the To, t, and 
P output from the fltting. Labeled on each plot is the IRDC and core name, classiflcation, 
and the derived parameters. 
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Fig. A-11. — 3 Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G025.04-00.20 (contour levels are 30, 60, 90, 120, 240, 360, 480, 840, 1200 mJy beam"^). 
The lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived 
from the millimeter, sub-millimeter, and far-lR continuum data are shown as filled circles 
(with the corresponding error bars), while the 24 fim fluxes are shown as either a fllled circle 
(when included within the fit), an open circle (when excluded from the fit), or as an up- 
per limit arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter 
regime (i.e. a limit at 24 /xm), we show the results from two flts: one using only the mea- 
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Fig. A-12. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G027.75+00.16 (contour levels are 30, 60, 90, 120, 180 mJy beam"^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-lR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /xm fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fiuxes 
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Fig. A-13. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G028.10-00.45, G028.04-00.46, G027.97-00.42, and G027.94-00.47 (contour levels are 
30, 60, 90, 120 mJy beam~^). The lower panels show the broadband SEDs for cores within 
this IRDC. The fluxes derived from the millimeter, sub-millimeter, and far-IR continuum 
data are shown as filled circles (with the corresponding error bars), while the 24 /xm fluxes 
are shown as either a filled circle (when included within the fit), an open circle (when ex- 
cluded from the fit), or as an upper limit arrow. For cores that have measured fluxes only in 
the millimeter/sub- millimeter regime (i.e. a limit at 24 /.im), we show the results from two 
fits: one using only the measured fluxes (solid line; lower limit), while the other includes the 
24 /im limit as a real data (dashed line; upper limit). In all other cases, the solid line is the 
best fit gray-body, while the dotted lines correspond to the functions determined using the 
errors for the Tjj, r, and /3 output from the fitting. Labeled on each plot is the IRDC and 
core name, classification, and the derived parameters. 
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Fig. A- 14. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G028. 08+00. 07 (contour levels are 30, 60, 90, 120 mJy beam~^). The lower panels show 
the broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, 
sub-millimeter, and far-IR continuum data are shown as filled circles (with the correspond- 
ing error bars), while the 24 /im fluxes are shown as either a filled circle (when included 
within the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For 
cores that have measured fluxes only in the millimeter /sub- millimeter regime (i.e. a limit 
at 24 /xm), we show the results from two fits: one using only the measured fluxes (solid line; 
lower limit), while the other includes the 24 /im limit as a real data (dashed line; upper 
limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
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Fig. A- 15. — Spitzer 24yum image overlaid with 1.2 mm continuum emission for 
G028. 23— 00.19 (contour levels are 30, 60, 90mJy beam^^). The lower panels show the 
broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub- 
millimeter, and far-IR continuum data are shown as filled circles (with the corresponding 
error bars), while the 24 /im fluxes are shown as either a filled circle (when included within 
the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For cores that 
have measured fluxes only in the millimeter/sub- millimeter regime (i.e. a limit at 24 /im), 
we show the results from two fits: one using only the measured fluxes (solid line; lower limit), 
while the other includes the 24 yum limit as a real data (dashed line; upper limit). In all 
other cases, the solid line is the best fit gray-body, while the dotted lines correspond to the 
functions determined using the errors for the T^, r, and /3 output from the fitting. Labeled 
on each plot is the IRDC and core name, classification, and the derived parameters. 
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Fig. A-16. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G028.28-00.34 (contour levels are 30, 60, 90, 120, 240, 360, 480 mJy beam"^). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /xm fluxes are shown as either a flUed circle (when 
included within the flt), an open circle (when excluded from the flt), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /im), we show the results from two flts: one using only the measured fluxes 
(solid line; lower limit), while the other includes the 24 //m limit as a real data (dashed line; 
upper limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the Td, t, and /3 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A- 17. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G028.37+00.07 (contour levels are 30, 60, 90, 120, 240, 360, 480, 840, 1200 mJy beam"^). 
The lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived 
from the millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles 
(with the corresponding error bars), while the 24 fim fiuxes are shown as either a filled circle 
(when included within the fit), an open circle (when excluded from the fit), or as an up- 
per limit arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter 
regime (i.e. a limit at 24 //m), we show the results from two fits: one using only the mea- 
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Fig. A- 18. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G028.53-00.25 (contour levels are 30, 60, 90,120, 240 mJy beam-^). The lower panels show 
the broadband SEDs for cores within this IRDC. The fluxes derived from the milhmeter, 
sub-milhmeter, and far-lR continuum data are shown as filled circles (with the correspond- 
ing error bars), while the 24 /xm fluxes arc shown as either a filled circle (when included 
within the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For 
cores that have measured fluxes only in the millimeter/sub- millimeter regime (i.e. a limit 
at 24 //m), we show the results from two flts: one using only the measured fluxes (solid line; 
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Fig. A- 19. — Spitzer 2Afj,m image overlaid with 1.2 mm continuum emission for 
G028. 67+00. 13 (contour levels are 30, 60, 90mJy beam^^). The lower panels show the 
broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub- 
millimeter, and far-IR continuum data are shown as filled circles (with the corresponding 
error bars), while the 24 /im fluxes are shown as either a flUed circle (when included within 
the flt), an open circle (when excluded from the flt), or as an upper limit arrow. For cores that 
have measured fluxes only in the millimeter/sub- millimeter regime (i.e. a limit at 24 /im), 
we show the results from two fits: one using only the measured fluxes (solid line; lower limit), 
while the other includes the 24 yum limit as a real data (dashed line; upper limit). In all 
other cases, the solid line is the best fit gray-body, while the dotted lines correspond to the 
functions determined using the errors for the T^, r, and /3 output from the fitting. Labeled 
on each plot is the IRDC and core name, classification, and the derived parameters. 
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Fig. A-20. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G030. 14—00.06 (contour levels are 30, 60mJy beam~^). The lower panels show the broad- 
band SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub- 
millimeter, and far-IR continuum data are shown as filled circles (with the corresponding 
error bars), while the 24 fim fluxes are shown as either a filled circle (when included within 
the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For cores that 
have measured fluxes only in the millimeter/sub- millimeter regime (i.e. a limit at 24yum), 
we show the results from two fits: one using only the measured fluxes (solid line; lower limit). 
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Fig. A-21. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G030.57-00.23 (contour levels are 30, 60, 90, 120, 240, 360, 480 mJy beam-^). The lower 
panels show the broadband SEDs for the cores within this IRDC. The lower panels show 



-77- 



-1°44' 



O 
o 
o 



Pi 

o 
Q 




Right Ascension (J2000) 




Fig. A-22. — Spitzer 24yum image overlaid with 1.2 mm continuum emission for 
G031.02-00.10 and G030.97-00.14 (contour levels are 40, 60, 80mJy beam'^). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /im fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter /sub-millimeter regime 
(i.e. a limit at 24 yUm), we show the results from two fits: one using only the measured fluxes 
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Fig. A-23. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G031.97+00.07 (contour levels are 30, 60, 90,120, 240, 360, 480 mJy beam-^). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24/i,m fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fluxes 
(solid line; lower limit), while the other includes the 24 //m limit as a real data (dashed line; 
upper limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the To, t, and (3 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A-24. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G033.69-00.01 (contour levels are 30, 60, 90, 120, 240 mJy beam-^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-lR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 /xm fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fiuxes 
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Fig. A-25. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G034.43+00.24 (contour levels are 60, 90, 120, 240, 360, 480, 840, 1200 mJy beam-^). The 
lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived from 
the millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with 
the corresponding error bars), while the 24 /xm fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fiuxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fiuxes 
(solid line; lower limit), while the other includes the 24 //m limit as a real data (dashed line; 
upper limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the To, t, and (3 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A-26. — Spitzer 24yum image overlaid with 1.2 mm continuum emission for 
G034. 77— 00.55 (contour levels are 30, 60, 90, 120, 240 mJy beam~^). The lower panels 
show the broadband SEDs for cores within this IRDC. The fluxes derived from the mil- 
limeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24yum fluxes are shown as either a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub- millimeter regime 
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Fig. A-27. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for G035.59- 
00.24 and G035.39-00.33 (contour levels are 30, 60, 90, 120 mJy beam"!). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24yum fiuxcs arc shown as cither a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fluxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 //m), we show the results from two fits: one using only the measured fluxes 
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Fig. A-28. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G036. 67— 00.11 (contour levels are 30, 60, 90, 120 mJy beam~^). The lower panels show 
the broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, 
sub-millimeter, and far-IR continuum data are shown as filled circles (with the correspond- 
ing error bars), while the 24 /im fluxes are shown as either a filled circle (when included 
within the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For 
cores that have measured fluxes only in the millimeter /sub- millimeter regime (i.e. a limit 
at 24 /im), we show the results from two fits: one using only the measured fluxes (solid line; 
lower limit), while the other includes the 24 /im limit as a real data (dashed line; upper 
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Fig. A-29. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G038.95— 0.47 (contour levels are 30, 60, 90, 120 mJy beam~^). The lower panels show 
the broadband SEDs for cores within this IRDC. The fluxes derived from the miUimeter, 
sub-millimeter, and far-IR continuum data are shown as filled circles (with the correspond- 
ing error bars), while the 24 //m fluxes are shown as either a filled circle (when included 
within the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For 
cores that have measured fiuxes only in the millimeter/sub- millimeter regime (i.e. a limit 
at 24 yum), we show the results from two fits: one using only the measured fiuxes (solid line; 
lower limit), while the other includes the 24 //m limit as a real data (dashed line; upper 
hmit). In all other cases, the solid fine is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the To, t, and (5 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A-30. — Spitzer 24 /xm image overlaid with 1.2 mm continuum emission for 
G048. 65— 00.29 (contour levels are 30, 60, 90mJy beam~^). The lower panels show the 
broadband SEDs for cores within this IRDC. The fluxes derived from the millimeter, sub- 
millimeter, and far-IR continuum data are shown as filled circles (with the corresponding 
error bars), while the 24 fim fluxes are shown as either a filled circle (when included within 
the fit), an open circle (when excluded from the fit), or as an upper limit arrow. For cores that 
have measured fluxes onlv in the millimeter /sub- millimeter reeime (i.e. a limit at 24umV 
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Fig. A-31. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G053.11+00.05 (contour levels are 30, 60, 90, 120, 240, 360, 480 mJy beam-^). The lower 
panels show the broadband SEDs for cores within this IRDC. The fluxes derived from the 
millimeter, sub-millimeter, and far-lR continuum data are shown as filled circles (with the 
corresponding error bars), while the 24 //m fluxes are shown as cither a filled circle (when 
included within the fit), an open circle (when excluded from the fit), or as an upper limit 
arrow. For cores that have measured fiuxes only in the millimeter/sub-millimeter regime 
(i.e. a limit at 24 /xm), we show the results from two fits: one using only the measured fiuxes 
(solid line; lower limit), while the other includes the 24 //m limit as a real data (dashed line; 
upper limit). In all other cases, the solid line is the best fit gray-body, while the dotted lines 
correspond to the functions determined using the errors for the To, t, and (5 output from 
the fitting. Labeled on each plot is the IRDC and core name, classification, and the derived 
parameters. 
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Fig. A-32. — Spitzer 24 //m image overlaid with 1.2 mm continuum emission for 
G053.31+00.00 and G053.25+00.04 (contour levels are 30, 60, 90, 120, 240 mJy beam-^). 
The lower panels show the broadband SEDs for cores within this IRDC. The fluxes derived 
from the millimeter, sub-millimeter, and far-IR continuum data are shown as filled circles 
(with the corresponding error bars), while the 24 ^m fluxes arc shown as either a fllled circle 
(when included within the flt), an open circle (when excluded from the fit), or as an up- 
per limit arrow. For cores that have measured fluxes only in the millimeter /sub- millimeter 
regime (i.e. a limit at 24 /xm), we show the results from two fits: one using only the mea- 
sured fiuxes (solid line; lower limit), while the other includes the 24 /im limit as a real data 
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Table A-1. Core coordinates, distances, and measured fluxes 
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Table A- 1— Continued 
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4.7 


< 0.06 






0.60 


1.18 


4.61 


1.03 


0.68 




0.06 


G028 


10-00 


45 


MMl 


18:44:12 


9 


-04:29:45 


3.2 


< 0.06 












0.25 






0.03 


G028 


10-00 


45 


MM2 


18:44:14 


3 


-04:29:48 


3.2 


< 0.05 












0.45 






0.02 


G028 


23-00 


19 


MMl 


18:43:30 


7 


-04:13:12 


4.8 


< 0.04 








4.69 


1.48 


1.67 


0.97 




0.08 


G028 


23-00 


19 


MM2 


18:43:29 





-04:12:16 


4.8 


< 0.06 












1.01 


0.30 




0.04 


G028 


23-00 


19 


MM3 


18:43:30 





-04:12:33 


4.8 


< 0.05 














0.20 




0.02 


G028 


28-00 


34 


MMl 


18:44:15 





-04:17:54 


3.2 


2.00 












15.66 


11.06 




0.32 


G028 


28-00 


34 


MM2 


18:44:21 


3 


-04:17:37 


3.2 


2.00 












12.00 


6.54 




0.17 


G028 


28-00 


34 


MM3 


18:44:13 


4 


-04:18:05 


3.2 


1.50 












7.67 


7.55 




0.14 


G028 


28-00 


34 


MM4 


18:44:11 


4 


-04:17:22 


3.2 


< 0.18 


















0.06 


G028 


37+00 


07 


MMl 


18:42:52 


1 


-03:59:45 


4.8 


0.12 












75.35 


37.33 




0.47 


G028 


37+00 


07 


MMIO 


18:42:54 





-04:02:30 


4.8 


0.02 












4.13 


2.59 




0.09 


G028 


37+00 


07 


MMll 


18:42:42 


7 


-04:01:44 


4.8 


0.01 












1.94 


1.21 




0.08 


G028 


37+00 


07 


MM12 


18:43:09 


9 


-04:06:52 


4.8 


< 0.07 


















0.07 


G028 


37+00 


07 


MM13 


18:42:41 


8 


-03:57:08 


4.8 


0.31 


















0.07 


G028 


37+00 


07 


MM14 


18:42:52 


6 


-04:02:44 


4.8 


< 0.04 












0.85 


1.00 




0.06 


G028 


37+00 


07 


MM15 


18:42:32 


4 


-04:01:16 


4.8 


0.28 


















0.06 


G028 


37+00 


07 


MM16 


18:42:40 


2 


-04:00:23 


4.8 


< 0.05 












1.81 


1.37 




0.06 


G028 


37+00 


07 


MM17 


18:43:00 





-04:01:34 


4.8 


< 0.04 












0.72 


0.36 




0.06 


G028 


37+00 


07 


MM18 


18:43:12 


9 


-04:01:16 


4.8 


0.43 


















0.06 


G028 


37+00 


07 


MM2 


18:42:37 


6 


-04:02:06 


4.8 


1.74 












26.54 


12.66 




0.28 


G028 


37+00 


07 


MM3 


18:43:03 


1 


-04:06:24 


4.8 


0.16 


















0.26 


G028 


37+00 


07 


MM4 


18:42:50 


7 


-04:03:16 


4.8 


0.04 












7.96 


6.48 




0.20 


G028 


37+00 


07 


MM5 


18:42:26 


8 


-04:01:30 


4.8 


< 0.00 


















0.15 


G02S 


37+00 


07 


MM6 


18:42:49 





-04:02:23 


4.8 


< 0.04 












5.22 


3.76 




0.15 


G028 


37+00 


07 


MM7 


18:42:56 


3 


-04:07:31 


4.8 


0.21 


















0.14 


G028 


37+00 


07 


MM8 


18:42:49 


7 


-04:09:64 


4.8 


1.06 


















0.12 


G028 


37+00 


07 


MM9 


18:42:46 


7 


-04:04:08 


4.8 


< 0.04 












2.24 


2.06 




0.10 


G028 


63-00 


25 


MMl 


18:44:18 





-03:59:34 


5.2 


< 0.06 








2.03 


3.88 


5.26 


4.90 


1.06 


0.23 


G028 


53-00 


25 


MMIO 


18:44:18 


5 


-03:58:43 


5.2 


< 0.06 














0.64 




0.06 


G028 


53-00 


25 


MM2 


18:44:15 


7 


-03:59:41 


5.2 


0.02 












2.37 


2.07 


0.48 


0.13 


G()28 


53-00 


25 


MM3 


18:44:16 





-04:00:48 


5.2 


< 0.05 




0.84 


1.66 


4.70 


11.03 


2.30 


1.55 


0.38 


0.13 


G028 


53-00 


25 


MM4 


18:44:18 


6 


-04:00:05 


5.2 


0.02 












1.50 


0.82 




0.10 


G028 


63-00 


25 


MM5 


18:44:17 





-04:02:04 


5.2 


0.01 














2.28 


0.67 


0.10 


G028 


53-00 


25 


MM6 


18:44:17 


8 


-04:00:06 


5.2 


< 0.06 












0.85 


0.87 




0.07 


G028 


53-00 


25 


MM7 


18:44:23 


7 


-04:02:09 


5.2 


< 0.06 






0.76 


1.08 


1.51 




0.80 


0.13 


0.06 


G028 


53-00 


25 


MM8 


18:44:22 





-04:01:35 


5.2 


< 0.06 












0.47 




0.20 


0.06 


G028 


53-00 


25 


MM9 


18:44:19 


3 


-03:58:05 


5.2 


< 0.05 












0.83 




0.29 


0.06 


G02S 


67+00 


13 


MMl 


18:43:03 


1 


-03:41:41 


4.8 


2.51 












8.66 


4.27 




0.09 


G028 


67+00 


13 


MM2 


18:43:07 


1 


-03:44:01 


4.8 


0.04 


0.64 


1.19 


2.63 


4.37 


7.63 


1.60 


1.37 




0.07 


G028 


67+00 


13 


MM3 


18:42:68 


2 


-03:48:20 


4.8 


0.36 


















0.06 


G028 


67+00 


13 


MM4 


18:43:13 


2 


-03:41:03 


4.8 


0.01 


















0.04 


G028 


67+00 


13 


MM5 


18:43:10 


1 


-03:45:08 


4.8 


< 0.05 












0.90 


0.39 




0.03 


G028 


67+00 


13 


MM6 


18:43:12 


2 


-03:45:39 


4.8 


< 0.06 












0.71 


0.33 




0.03 


G028 


67+00 


13 


MM7 


18:43:06 


9 


-03:45:11 


4.8 


< 0.05 


















0.03 


G03() 


14-00 


06 


MMl 


18:46:35 


7 


-02:31:03 


5.2 


< 0.07 


1.07 


1.62 


4.07 


6.89 


16.87 


0.95 


0.77 




0.05 


G030 


14-00 


06 


MM2 


18:46:31 


7 


-02:32:41 


5.2 


< 0.01 












1.07 


0.52 




0.04 


G030 


67-00 


23 


MMl 


18:48:00 





-02:07:20 


5.2 


0.67 


















0.30 


G030 


57-00 


23 


MM2 


18:47:58 


7 


-02:15:20 


5.2 


1.31 


















0.14 


G030 


57-00 


23 


MM3 


18:47:54 


5 


-02:11:15 


5.2 


0.01 


2.41 


2.91 


6.79 


10.67 


15.24 








0.06 


G03() 


57-00 


23 


MM4 


18:48:01 


8 


-02:12:35 


5.2 


< 0.06 


1.38 


2.27 


5.17 


5.15 


4.57 


1.25 






0.04 


G03() 


57-00 


23 


MM5 


18:47:59 


9 


-02:11:01 


5.2 


< 0.08 


















0.04 


G03() 


57-00 


23 


MM6 


18:47:59 


4 


-02:13:13 


5.2 


0.01 


















0.03 


G030 


97-00 


14 


MMl 


18:48:21 


6 


-01:48:27 


4.8 


1.66 


106.60 


127.42 


136.46 


130.72 


122.59 


16.44 


11.21 


1.23 


0.11 


G030 


97-00 


14 


MM2 


18:48:22 





-01:47:42 


4.8 


< 0.07 


















0.03 


G031 


02-00 


10 


MMl 


18:48:09 


9 


-01:45:17 


4.8 


< 0.08 


















0.13 


G031 


97+00 


07 


MMl 


18:49:36 


3 


00:45:45 


6.1 


0.86 












48.79 


40.60 


5.49 


0.91 


G031 


97+00 


07 


MM2 


18:49:36 





00:46:16 


6.1 


< 0.07 












9.78 


10.68 


1.08 


0.31 


G031 


97+00 


07 


MM3 


18:49:32 


3 


00:47:02 


6.1 


< 0.07 












5.67 


4.70 


0.75 


0.19 


G031 


97+00 


07 


MMl 


18:49:33 





00:47:33 


6.1 


0.02 












2.18 


2.45 


0.41 


0.12 


G031 


97+00 


07 


MM5 


18:49:21 


9 


00:50:35 


6.1 


0.01 












1.89 


2.35 


0.47 


0.10 


G031 


97+00 


07 


MM6 


18:49:35 





00:46:44 


6.1 


< 0.08 












2.97 


2.67 


0.66 


0.10 


G031 


97+00 


07 


MM7 


18:49:28 


4 


00:48:64 


6.1 


< 0.04 












0.87 


1.08 


0.46 


0.08 


G031 


97+00 


07 


MM8 


18:49:29 


1 


00:48:12 


6.1 


0.21 












1.44 


1.07 




0.06 


G031 


97+00 


07 


MM9 


18:49:31 


6 


00:46:30 


6.1 


< 0.07 


















0.06 
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Tabic A-1 — Contimicd 



IRDC core 


Coordinates 


D 


F24 


F55 


F65 


F75 


F85 


F95 


F350 


F450 


Fs50 


F1.2 




(J2000) (J2000) 


(kpc) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 


(Jy) 



G033 


69—00 


01 


MMl 


18 


52:58 


8 


00 


42 


37 


7.1 







21 


7.25 


12.68 


19.62 


25.87 


36.12 


14 


61 


13 


32 






0.21 


G033 


69—00 


01 


MMIO 


18 


52:52 


7 


00 


38 


35 


7.1 


< 





05 












1 


10 


1 


22 






0.04 


G033 


69—00 


01 


MMll 


18 


52:56 


2 


00 


41 


48 


7.1 


< 





05 












1 


10 


1 


17 






0.03 


G033 


69-00 


01 


MM2 


18 


52:49 


9 


00 


37 


57 


7.1 




2 


54 












8 


67 


6 


15 






0.12 


G033 


69-00 


01 


MM3 


18 


52:50 


8 


00 


36 


43 


7.1 







76 












3 


36 


2 


85 






0.08 


G033 


69-00 


01 


MM4 


18 


52:56 


4 


00 


43 


08 


7.1 







04 


1.06 


1.45 


2.38 


2.64 


4.67 


3 


70 


3 


68 






0.08 


G033 


69—00 


01 


MM5 


18 


52:47 


8 


00 


36 


47 


7.1 







97 












6 


07 


4 


52 






0.06 


G033 


69—00 


01 


1VIM6 


18 


52:48 


7 


00 


35 


58 


7.1 


< 





05 












1 


20 










t).t)5 


G033 


69-00 


01 


MM7 


18 


52:58 


1 


00 


44 


08 


7.1 


< 





05 
























0.04 


G033 


69-00 


01 


MM8 


18 


52:53 


9 


00 


41 


16 


7.1 


< 





05 












1 


36 


1 


04 






0.04 


G033 


69-00 


01 


MM9 


18 


52:58 


1 


00 


41 


20 


7.1 







57 












2 


70 


2 


08 






0.04 


G034 


43+00 


24 


MMl 


18 


53:18 





01 


25 


24 


3.7 




1 


46 












147 


14 


50 


85 


10 


10 


2.23 


G034 


43+00 


24 


MM 2 


18 


53:18 


6 


01 


24 


40 


3.7 




1 


62 












55 


09 


24 


98 


5 


51 


0.96 


G034 


43+00 


24 


MM3 


18 


53:20 


4 


01 


28 


23 


3.7 







12 


18.85 


31.21 


42.72 


5 1.51 


55.65 


12 


36 


9 


63 


1 


57 


0.24 


G034 


43+00 


24 


MM4 


18 


53:19 





01 


24 


08 


3.7 







13 


21.93 


44.68 


75.76 


100.36 


113.95 


5 


77 


4 


15 





84 


0.22 


G034 


43+00 


24 


MM5 


18 


53:19 


8 


01 


23 


30 


3.7 


< 





05 


0.86 


3.33 


6.27 


11.58 


16.85 


3 


13 


1 


41 





25 


0.12 


G034 


43+00 


24 


MM6 


18 


53:18 


6 


01 


27 


48 


3.7 


< 





04 















51 





56 





15 


0.06 


G034 


43+00 


24 


MM7 


18 


53:18 


3 


01 


27 


13 


3.7 


< 





05 












1 


52 





59 





20 


0.06 


G034 


43+00 


24 


MM8 


18 


53:16 


4 


01 


26 


20 


3.7 


< 





05 


4.23 


3.72 


2.34 


9.55 


21.63 





99 


1 


28 





16 


0.05 


G034 


43+00 


24 


MM9 


18 


53:18 


4 


01 


28 


14 


3.7 


< 





04 
























0.05 


G034 


77-00 


55 


MMl 


18 


56:48 


2 


01 


18 


47 


2.9 




2 


04 
























0.23 


G034 


77-00 


55 


MM2 


18 


56:50 


3 


01 


23 


16 


2.9 


< 





04 


1.79 


1.64 


2.62 


3.19 


6.30 














0.06 


G034 


77-00 


55 


MM3 


18 


56:44 


7 


01 


20 


42 


2.9 







00 
























0.04 


G034 


77-00 


55 


MM4 


18 


56:48 


9 


01 


23 


34 


2.9 


< 





04 
























0.04 


G035 


39-00 


33 


MMl 


18 


56:41 


2 


02 


09 


52 


3.0 


< 





00 
























0.13 


G035 


39-00 


33 


MM2 


18 


56:59 


2 


02 


04 


53 


3.0 







27 
























0.11 


G035 


39-00 


33 


MM3 


18 


57:05 


3 


02 


06 


29 


3.0 




2 


37 












11 


45 


6 


72 






0.11 


G035 


39-00 


33 


MM4 


18 


57:06 


7 


02 


08 


23 


3.0 







01 












1 


80 


1 


08 






0.08 


G035 


39-00 


33 


MM5 


18 


57:08 


8 


02 


08 


09 


3.0 


< 





04 


0.15 


0.81 


0.69 


2.16 


4.29 














0.07 


G035 


39-00 


33 


MM6 


18 


57:08 


4 


02 


09 


09 


3.0 







03 












2 


56 


1 


57 






0.06 


G035 


39-00 


33 


MM7 


18 


57:08 


1 


02 


10 


50 


3.0 







02 












2 


18 


1 


48 






0.06 


G035 


39-00 


33 


JvlJVlo 


18 


57:07 





02 


08 


54 


3.0 


< 





04 















33 










0.05 


G035 


39-00 


33 


MM9 


18 


57:11 


2 


02 


07 


27 


3.0 


< 





03 



















86 






0.05 


G035 


59-00 


24 


MMl 


18 


57:02 


3 


02 


17 


04 


3.0 







36 












6 


11 










0.07 


G035 


59-00 


24 


MM2 


18 


57:07 


4 


02 


16 


14 


3.0 







02 


2.01 


4.11 


7.49 


12.73 


17.87 


2 


17 










0.06 


G035 


59-00 


24 


MM3 


18 


57:11 


6 


02 


16 


08 


3.0 


< 





03 















65 










0.04 


G036 


67-00 


11 


MMl 


18 


58:39 


6 


03 


16 


16 


3.6 


< 





03 


0.11 


0.90 


1.59 


4.05 


5.14 


1 


52 










0.06 


G036 


67-00 


11 


MM2 


18 


58:35 


6 


03 


15 


06 


3.6 


< 





03 


0.29 


0.39 


1.71 


4.21 


9.51 














0.05 


G038 


95-00 


47 


MMl 


19 


04:07 


4 


05 


08 


48 


2.8 


< 





03 


0.02 


0.49 


0.55 


0.80 


3.84 


4 


04 


2 


81 





79 


0.12 


G038 


95-00 


47 


MM2 


19 


04:03 


4 


05 


07 


56 


2.8 




1 


25 












6 


97 


3 


53 





43 


0.10 


G038 


95-00 


47 


MM3 


19 


04:07 


4 


05 


09 


44 


2.8 







00 












3 


31 


1 


94 





33 


0.07 


G038 


95-00 


47 


MM4 


19 


04:00 


6 


05 


09 


06 


2.8 







26 












1 


81 


1 


24 





27 


0.05 


G048 


65-00 


29 


MMl 


19 


21:49 


7 


13 


49 


30 


2.6 







10 












4 


02 


2 


39 






0.07 


G048 


65-00 


29 


MM2 


19 


21:47 


6 


13 


49 


22 


2.6 


< 





04 












1 


19 


1 


34 






0.06 


G053 


11+00 


05 


MMl 


19 


29:17 


2 


17 


56 


21 


1.8 




1 


37 












49 


71 


22 


42 






0.54 


G053 


11+00 


05 


MM2 


19 


29:20 


2 


17 


57 


06 


1.8 







18 












4 


78 


2 


45 






0.08 


G053 


11+00 


05 


MM3 


19 


29:00 


6 


17 


55 


11 


1.8 







00 
























0.07 


G053 


11 + 00 


05 


MM4 


19 


29:20 


4 


17 


55 


04 


1.8 







23 












2 


76 


1 


09 






0.05 


G053 


11+00 


05 


MM5 


19 


29:26 


3 


17 


54 


53 


1.8 







63 












2 


11 










0.05 


G053 


25+00 


04 


MMl 


19 


29:39 





18 


01 


42 


2.0 







06 












8 


15 


4 


30 






0.25 


G053 


25+00 


04 


MM2 


19 


29:33 





18 


01 


00 


2.0 







07 












1 


25 





58 






0.10 


G053 


25+00 


04 


MM3 


19 


29:44 





17 


58 


47 


2.0 







01 












2 


70 





93 






0.10 


G053 


25+00 


04 


MM4 


19 


29:34 


5 


18 


01 


39 


2.0 







09 












3 


00 


1 


68 






0.07 


G053 


25+00 


04 


MM5 


19 


29:39 


4 


17 


58 


40 


2.0 







01 












1 


79 





65 






0.05 


G053 


25 + 00 


04 


MM6 


19 


29:31 




17 


59 


50 


2.0 




1 


39 












2 


58 


1 


13 






0.05 


G053 


31+00 


00 


MMl 


19 


29:50 





18 


05 


07 


2.0 


< 





03 


0.20 


0.38 


0.70 


0.77 


2.40 





73 





52 






0.07 


G053 


31+00 


00 


MM2 


19 


29:42 


1 


18 


03 


57 


2.0 


< 





03 


0.35 


1.28 


2.51 


3.30 


5.35 





78 





78 






0.06 


G053 


31+00 


00 


MM3 


19 


29:49 


7 


18 


04 


39 


2.0 







00 


0.95 


1.63 


2.51 


3.28 


5.18 


1 


58 





79 






0.02 
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Table A-2. Properties of the cores derived from gray-body fits to their broadband SEDs. 



irtuo core 


T' 
i^D 


''"260 


a 
P 


T 

Li 


IVl 






X 


Classification^ 














(cm ) 


(cm ) 






G015. 31-00. 16 MMl 


29 


0.004 


1.6 


56 


7 


22.41 


5.69 


1.91 


R 


G015. 31-00. 16 MM2 


23 


0.003 


1.5 


16 


21 


22.88 


6.16 


9.72 


I 


G018. 82-00. 28 MMl 


26 


0.091 


2.0 


4526 


495 


23.51 


6.42 


1.44 


I 


G018.82-00.28 MM2 


38 


0.020 


1.9 


8990 


87 


22.76 


5.66 


1.65 


R 


G018.82-00.28 MM3 


34 


0.034 


2.4 


11356 


103 


22.83 


5.73 


0.61 


A 


G018.82-00.28 MM4 


17 


0.031 


2.1 


198 


228 


23.18 


6.08 


4.94 


I 


G019. 27+00. 07 MMl 


33 


0.027 


2.1 


1998 


52 


23.09 


6.27 


1.50 


B 


G019.27+00.07 MM2 


30 


0.017 


2.1 


760 


60 


23.15 


6.33 


1.27 


A 


G022. 35+00. 41 MMl 


20 


0.028 


1.1 


138 


182 


23.27 


6.26 


9.37 


A 


G022. 35+00. 41 MM2 


44 


0.001 


0.9 


267 


60 


22.78 


5.78 


0.85 


R 


G022.73+00.11 MMl 


29 


0.001 


1.0 


80 


64 


22.57 


5.44 


5.13 


Q 


G023.60+00.00 MMl 


28 


0.024 


1.3 


764 


184 


23.27 


6.27 


4.72 


A 


G023. 60+00. 00 MM2 


30 


0.028 


1.6 


1788 


108 


23.04 


6.04 


2.04 


A 


G023. 60+00. 00 MM3 


46 


0.001 


1.0 


470 


79 


22.90 


5.90 


4.28 


R 


G023.60+00.00 MM5 


48 


0.002 


1.3 


1272 


27 


22.43 


5.43 


0.19 


A 


G023.60+00.00 MM7 


44 


0.001 


0.7 


137 


80 


22.91 


5.90 


1.56 


I 


G024.08+00.04 MMl 


37 


0.052 


2.4 


14847 


71 


22.86 


5.85 


0.02 


R 


G024.08+00.04 MM2 


29 


0.006 


1.6 


321 


115 


23.07 


6.06 


1.34 


Q 


G024.08+00.04 MM4 


30 


0.002 


1.0 


60 


55 


22.74 


5.74 


2.06 


Q 


G024. 33+00. 11 MMl 


37 


0.105 


1.7 


47387 


472 


23.21 


5.98 


0.90 


R 


G024.33+00.11 MMIO 


37 


0.002 


1.0 


187 


36 


22.55 


5.55 


< 0.01 


I 


G024.33+00.11 MM9 


52 


0.001 


1.0 


697 


81 


22.91 


5.91 


< 0.01 


R 


G024.60+00.08 MMl 


25 


0.038 


1.7 


964 


163 


23.23 


6.24 


2.84 


A 


G024.60+00.08 MM2 


23 


0.073 


2.1 


3544 


249 


22.91 


5.67 


1.21 


I 


G025.04-00.20 MMl 


34 


0.006 


1.0 


498 


189 


23.18 


6.13 


0.44 


A 


G025. 04-00.20 MM2 


28 


0.002 


0.5 


37 


95 


22.88 


5.83 


7.29 


I 


G027.75+00.16 MMl 


20 


0.015 


1.8 


224 


205 


23.09 


5.98 


1.64 


A 


G027. 75+00. 16 MM2 


24 


0.010 


2.0 


376 


80 


22.68 


5.56 


0.36 


Q 


G027.75+00.16 MM3 


40 


0.001 


0.7 


109 


86 


22.71 


5.59 


< 0.01 


I 


G027.75+00.16 MM5 


19 


0.009 


2.3 


108 


50 


22.48 


5.36 


4.17 


Q 


G027.94-00.47 MMl 


44 


0.004 


1.4 


1260 


23 


22.48 


5.54 


0.05 


R 


G027.97-00.42 MMl 


33 


0.012 


1.6 


890 


25 


22.52 


5.59 


1.19 


A 


G027.97-00.42 MM2 


34 


0.018 


2.1 


2743 


26 


22.54 


5.61 


0.40 


R 


G028. 04-00. 46 MMl 


38 


0.007 


1.6 


1160 


31 


22.61 


5.67 


0.96 


A 


G028.08+00.07 MMl 


22 


0.004 


1.1 


47 


237 


23.18 


6.08 


6.02 


Q 


G028.23-00.19 MMl 


22 


0.004 


1.0 


61 


567 


23.53 


6.41 


6.37 


Q 


G028.28-00.34 MMl 


40 


0.020 


1.6 


4768 


150 


23.29 


6.35 


0.53 


R 


G028.28-00.34 MM2 


40 


0.015 


1.8 


4565 


132 


23.24 


6.29 


0.01 


R 


G028.28-00.34 MM3 


40 


0.010 


1.7 


2951 


23 


22.48 


5.53 


1.53 


R 


G028.37+00.07 MMl 


33 


0.165 


2.5 


26062 


507 


23.49 


6.37 


0.40 


A 


G028.37+00.07 MMIO 


32 


0.007 


1.6 


1200 


172 


23.02 


5.91 


0.22 


A 


G028.37+00.07 MMll 


36 


0.002 


1.0 


353 


173 


23.02 


5.91 


0.05 


A 


G028.37+00.07 MM2 


39 


0.041 


2.1 


23329 


184 


23.05 


5.93 


0.07 


I 


G028.37+00.07 MM4 


32 


0.014 


1.5 


2113 


140 


22.93 


5.82 


0.74 


A 


G028.53-00.25 MMl 


17 


0.038 


1.5 


211 


1088 


23.75 


6.60 


1.53 


Q 
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Table A- 2 — Continued 



IRDC core 


Tn 


T250 


B 


L 


M 






A 






CK') 


















G028. 53-00. 25 MM2 


38 


0.003 


0.8 


502 


709 


23.56 


6.41 


1.11 


A 


G028. 53-00. 25 MM3 


22 


0.009 


1.2 


163 


1225 


23.80 


6.65 


6.39 


Q 


G028. 53-00. 25 MM4 


41 


0.001 


0.5 


229 


270 


23.14 


5.99 


0.02 


I 


G028.53-00.25 MM5 


30 


0.010 


1.7 


1283 


104 


22.73 


5.58 


1.20 


I 


G028.53-00.25 MM7 


23 


0.003 


1.0 


58 


359 


23.27 


6.11 


0.24 


Q 


G028.67+00.13 MMl 


40 


0.013 


2.1 


11997 


46 


22.43 


5.31 


0.12 


R 


G028.67+00.13 MM2 


23 


0.006 


1.4 


145 


245 


23.16 


6.05 


2.64 


I 


G030. 14-00.06 MMl 


27 


0.003 


1.3 


163 


294 


23.18 


6.03 


9.35 


I 


G030. 57-00. 23 MM3 


23 


0.011 


1.9 


357 


161 


22.92 


5.77 


3.34 


I 


G030. 57-00. 23 MM4 


29 


0.002 


1.2 


162 


164 


22.92 


5.77 


1.46 


Q 


G030.97-00.14 MMl 


37 


0.012 


1.8 


4416 


179 


23.03 


5.92 


7.76 


A 


G031.97+00.07 MMl 


37 


0.081 


1.7 


36427 


574 


23.34 


6.12 


2.12 


A 


G031.97+00.07 MM4 


37 


0.003 


0.8 


638 


280 


23.02 


5.81 


1.60 


I 


G031. 97+00.07 MM5 


35 


0.003 


0.9 


534 


57 


22.33 


5.11 


2.73 


I 


G031. 97+00. 07 MM8 


47 


0.001 


0.9 


871 


124 


22.67 


5.45 


0.26 


A 


G033. 69-00. 01 MMl 


25 


0.023 


1.6 


2003 


751 


23.32 


6.03 


4.61 


R 


G033. 69-00. 01 MM2 


41 


0.013 


1.9 


22373 


563 


23.19 


5.90 


0.75 


R 


G033.69-00.01 MM3 


43 


0.004 


1.4 


5108 


98 


22.43 


5.15 


0.90 


R 


G033.69-00.01 MM4 


28 


0.003 


0.9 


246 


485 


23.13 


5.84 


7.56 


A 


G033.69-00.01 MM5 


37 


0.011 


2.3 


18213 


99 


22.44 


5.15 


1.08 


A 


G033.69-00.01 MM9 


39 


0.004 


1.8 


5621 


43 


22.07 


4.79 


0.87 


I 


G034.43+00.24 MMl 


38 


0.169 


1.6 


24360 


433 


23.64 


6.64 


0.74 


A 


G034.43+00.24 MM2 


38 


0.071 


1.6 


14304 


533 


23.73 


6.73 


1.10 


R 


G034.43+00.24 MM3 


29 


0.020 


1.5 


934 


171 


23.24 


6.24 


3.47 


A 


G034.43+00.24 MM4 


32 


0.015 


1.5 


1164 


115 


23.06 


6.06 


6.09 


A 


G034.43+00.24 MM5 


23 


0.013 


1.6 


177 


388 


23.60 


6.59 


5.76 


I 


G034.43+00.24 MM8 


43 


0.001 


0.5 


79 


39 


22.60 


5.59 


19.81 


I 


G034.77-00.55 MM2 


34 


0.001 


0.7 


40 


84 


23.13 


6.22 


3.75 


Q 


G035. 39-00. 33 MM3 


39 


0.019 


2.2 


6495 


32 


22.68 


5.77 


0.46 


R 


G035.39-00.33 MM4 


35 


0.002 


1.0 


116 


52 


22.90 


5.99 


0.02 


A 


G035.39-00.33 MM5 


18 


0.017 


1.8 


45 


129 


23.29 


6.38 


5.77 


Q 


G035.39-00.33 MM6 


34 


0.004 


1.6 


327 


47 


22.85 


5.94 


0.18 


A 


G035.39-00.33 MM7 


34 


0.004 


1.5 


242 


52 


22.90 


5.99 


0.29 


A 


G035. 59-00. 24 MMl 


36 


0.010 


2.0 


1685 


31 


22.67 


5.76 


< 0.01 


R 


G035. 59-00. 24 MM2 


24 


0.008 


1.8 


128 


31 


22.67 


5.76 


1.83 


A 


G036.67-00.il MMl 


20 


0.010 


1.7 


69 


37 


22.59 


5.60 


4.44 


Q 


G036.67-00.il MM2 


16 


0.047 


2.5 


137 


71 


22.88 


5.89 


4.81 


Q 


G038.95-00.47 MMl 


16 


0.033 


1.7 


40 


156 


23.43 


6.54 


9.49 


I 


G038.95-00.47 MM2 


40 


0.009 


1.8 


2031 


30 


22.70 


5.82 


0.20 


R 


G038.95-00.47 MM3 


28 


0.008 


1.7 


240 


7 


22.04 


5.15 


0.71 


I 


G038. 95-00.47 MM4 


43 


0.002 


1.2 


320 


25 


22.63 


5.74 


1.09 


A 


G048.65-00.29 MMl 


35 


0.007 


1.7 


537 


26 


22.73 


5.88 


0.21 


A 


G053.11+00.05 MMl 


38 


0.076 


2.1 


4473 


45 


23.28 


6.59 


0.02 


R 


G053.11+00.05 MM2 


36 


0.007 


1.7 


325 


19 


22.90 


6.21 


0.04 


A 


G053.11+00.05 MM4 


40 


0.003 


1.5 


178 


19 


22.89 


6.20 


0.10 


A 
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Table A- 2 — Continued 



IRDC core 


To 


1"2B0 





L 


M 


Log[N(H2)] 


Log[n(H2)] 




Classification'' 




(K) 








(Mo) 


(cm-2) 


(cm-3) 






G053.11+00.05 MM5 


45 


0.002 


1.4 


192 


5 


22.26 


5.57 


< 0.01 


R 


G053. 25+00. 04 MMl 


35 


0.010 


1.2 


306 


11 


22.55 


5.81 


< 0.01 


I 


G053. 25+00. 04 MM2 


47 


0.001 


0.5 


46 


24 


22.91 


6.18 


0.77 


B 


G053.25+00.04 MM3 


35 


0.002 


0.9 


49 


5 


22.16 


5.42 


0.65 


I 


G053.25+00.04 MM4 


37 


0.004 


1.4 


188 


6 


22.25 


5.51 


0.06 


A 


G053.25+00.04 MM5 


35 


0.002 


1.1 


46 


12 


22.58 


5.85 


0.36 


I 


G053.25+00.04 MM6 


46 


0.002 


1.5 


393 


7 


22.33 


5.59 


< 0.01 


A 


G053.31+00.00 MMl 


28 


0.001 


0.5 


7 


20 


22.82 


6.09 


3.39 


Q 


G053. 31+00. 00 MM2 


25 


0.003 


1.2 


17 


17 


22.77 


6.03 


5.67 


Q 


G053.31+00.00 MM3 


25 


0.003 


1.6 


21 


8 


22.42 


5.68 


2.12 


I 



denotes the classification based on the IRAC colors and the presence or absence of 24 //m emission. 
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Table A-3. Properties of the cores for which hmits where estimated from the gray-body 

fits to their broadband SEDs. 



IRDC core 


(K) 


''"260 






L 

{Lq 


■) 


M 






Classification'^ 


G015. 05+00. 07 


MMl 


11 - 


36 


0.190 - 


0.006 


2.5 - 


1.4 


22 - 


515 


225 - 


50 


0.68 - 


1.27 


O 


G015. 05+00. 07 


MM2 


10 - 


51 


0.012 - 


0.000 


1.2 - 


0.5 


1 


- 57 


209 - 


27 


< 0.01 


- 2.45 


R 


G015. 05+00. 07 


MM4 


14 - 


47 


0.010 - 


0.001 


1.6 - 


0.8 


g 


- 96 


61 - 


15 


< 0.01 - 


< 0.01 


O 


G015. 05+00. 07 


MM5 


13 - 


49 


0.010 - 


0.000 


1.6 - 


0.6 


4 


- 64 


47 - 


10 


< 0.01 - 


< 0.01 


o 


G015. 31-00. 16 


MM4 


13 - 


48 


0.010 - 


0.000 


1.8 - 


0.8 


2 


- 27 


13 - 


3 


< 0.01 - 


< 0.01 


o 


G015.31-00.16 


MM5 


15 - 


41 


0.010 - 


0.001 


2.1 - 


1.4 


4 


- 60 


24 - 


8 


< 0.01 - 


< 0.01 


o 


G018.82-00.28 


MM5 


18 - 


41 


0.010 - 


0.002 


1.6 - 


1.1 


42 - 


432 


169 - 


67 


0.21 - 


0.09 


o 


G022.35+00.41 


MM3 


19 - 


38 


0.010 - 


0.003 


1.6 - 


1.4 


38 - 


457 


101 - 


48 


< 0.01 - 


< 0.01 


B 


G023. 60+00. 00 


MM4 


11 - 


41 


0.114 - 


0.003 


2.5 - 


1.3 


17 - 


629 


729 - 


135 


0.02 - 


0.26 


R 


G023. 60+00. 00 


MM6 


12 - 


51 


0.011 - 


0.001 


1.2 - 


0.5 


3 _ 


125 


181 - 


30 


< 0.01 


- 0.26 


O 


G023. 60+00. 00 


MM8 


13 - 


54 


0.010 - 


0.000 


1.5 - 


0.5 


5 _ 


112 


65 - 


13 


< 0.01 - 


< 0.01 


O 


G023. 60+00. 00 


MM9 


13 - 


52 


0.010 - 


0.000 


1.5 - 


0.5 


4 


- 95 


135 - 


25 


0.13 - 


0.02 


o 


G024.08+00.04 


MM3 


20 - 


38 


0.010 - 


0.003 


1.8 - 


1.4 


62 - 


427 


101 - 


53 


0.01 - < 0.01 


o 


G024.08+00.04 


MM5 


18 - 


40 


0.010 - 


0.002 


1.8 - 


1.3 


30 - 


311 


89 - 


37 


< 0.01 


- 0.01 


o 


G024.33+00.11 


MMll 


9 - 


43 


0.108 - 


0.001 


2.5 - 


0.9 


6 - 


213 


975 - 


129 


0.38 - 


1.01 


o 


G024. 33+00.11 


MM2 


19 - 


46 


0.010 - 


0.001 


1.2 - 


0.6 


32 - 


269 


124 - 


47 


< 0.01 - 


< 0.01 


o 


G024. 33+00. 11 


MM3 


23 - 


40 


0.010 - 


0.003 


1.4 - 


1.2 


109 - 


526 


90 - 


53 


< 0.01 - 


< 0.01 


R 


G024. 33+00. 11 


MM4 


19 - 


43 


0.010 - 


0.002 


1.3 - 


0.8 


35 - 


286 


427 - 


174 


< 0.01 - 


< 0.01 


O 


G024.33+00.11 


MM5 


18 - 


46 


0.010 - 


0.001 


1.4 - 


0.8 


26 - 


284 


198 - 


70 


< 0.01 - 


< 0.01 


o 


G024.33+00.11 


MM6 


21 - 


40 


0.010 - 


0.002 


1.5 - 


1.1 


67- 


410 


230 - 


119 


< 0.01 - 


< 0.01 


o 


G024. 33+00.11 


MM8 


16 - 


47 


0.010 - 


0.001 


1.3 - 


0.6 


15 - 


182 


369 - 


108 


< 0.01 - 


< 0.01 


o 


G024. 60+00. 08 


MM3 


27 - 


42 


0.010 - 


0.004 


0.9 - 


0.7 


164 - 


511 


129 - 


87 


< 0.01 - 


< 0.01 


I 


G025. 04-00. 20 


MM3 


10 - 


41 


0.198 - 


0.002 


2.5 - 


1.0 


14 - 


420 


191 - 


28 


0.83 - 


1.62 


o 


G025.04-00.20 


MM4 


19 - 


40 


0.010 - 


0.002 


1.4 - 


1.0 


48 - 


405 


439 - 


197 


< 0.01 - 


< 0.01 


o 


G025.04-00.20 


MM5 


8 - 


48 


0.233 - 


0.001 


2.5 - 


0.5 


5 - 


164 


317 - 


27 


0.74- 


1.79 


o 


G027.75+00.16 


MM4 


20 - 


36 


0.010 - 


0.003 


2.1 - 


1.7 


123 - 


794 


95 - 


53 


< 0.01 - 


< 0.01 


o 


G027.94-00.47 


MM2 


10 - 


52 


0.010 - 


0.000 


1.2 - 


0.5 


1 


- 49 


127 - 


15 


< 0.01 


- 2.19 


o 


G028. 10-00.45 


MMl 


10 - 


51 


0.010 - 


0.000 


1.3 - 


0.5 


1 


- 45 


199 - 


26 


< 0.01 


- 0.88 


o 


G028. 10-00.45 


MM2 


12 - 


46 


0.010 - 


0.000 


1.9 - 


0.9 


3 


- 65 


40 - 


8 


< 0.01 - 


< 0.01 





G028. 23-00.19 


MM2 


14 - 


45 


0.010 - 


0.001 


1.7 - 


0.8 


11 - 


221 


139 - 


34 


0.94 - 


0.58 


I 


G028.23-00.19 


MM3 


10 - 


52 


0.010 - 


0.000 


1.6 - 


0.5 


3 


- 93 


87- 


12 


< 0.01 - 


< 0.01 


Q 


G028.37+00.07 


MM14 


8 - 


47 


0.185 - 


0.001 


2.5 - 


0.5 


7- 


190 


123 - 


12 


0.25 - 


0.98 


Q 


G028.37+00.07 


MM16 


10 - 


39 


0.106 - 


0.002 


2.5 - 


1.2 


20 - 


556 


1013 - 


185 


0.06 - 


0.41 


Q 


G028.37+00.07 


MM17 


13 - 


48 


0.010 - 


0.000 


1.2 - 


0.5 


7- 


157 


317 - 


63 


0.28- 


0.36 


Q 


G028. 37+00.07 


MM6 


11 - 


35 


0.235 - 


0.008 


2.5 - 


1.4 


60 - 1394 


449 - 


104 


0.06 - 


0.39 


I 


G028. 37+00.07 


MM9 


10 - 


39 


0.211 - 


0.003 


2.5 - 


1.0 


20 - 


546 


1173 - 


177 


0.17 - 


0.74 


Q 


G028. 53-00.25 


MMIO 


15 - 


45 


0.010 - 


0.001 


1.4 - 


0.7 


17 - 


291 


563 - 


152 


< 0.01 - 


< 0.01 


Q 


G028.53-00.25 


MM6 


8 - 


47 


0.231 - 


0.001 


2.5 - 


0.5 


7 - 


239 


413 - 


34 


0.03 - 


0.66 


Q 


G028.53-00.25 


MM8 


7- 


49 


0.262 - 


0.000 


2.5 - 


0.5 


4 - 


191 


989 - 


60 


0.07 - 


2.55 


Q 


G028. 53-00.25 


MM9 


8 - 


48 


0.188 - 


0.001 


2.5 - 


0.5 


7 - 


218 


662 - 


56 


0.66 - 


1.51 


Q 


G028. 67+00.13 


MM5 


14 - 


44 


0.010 - 


0.001 


1.7 - 


1.0 


14 - 


250 


238 - 


64 


0.23 - 


0.08 


Q 


G028. 67+00.13 


MM6 


13 - 


47 


0.010 - 


0.000 


1.7 - 


0.8 


9 - 


196 


173 - 


38 


0.18 - 


0.04 


Q 


G030.14-00.06 


MM2 


15 - 


37 


0.010 - 


0.001 


1.7- 


1.1 


24 - 


236 


99 - 


37 


0.09 - 


0.02 


Q 


G031.97+00.07 


MM2 


11 - 


35 


0.515 - 


0.016 


2.5 - 


1.3 


170 - 3970 


1566 - 


358 


0.81 - 


1.60 


Q 


G031.97+00.07 


MM3 


11 - 


36 


0.344 - 


0.008 


2.5 - 


1.2 


92 - 2443 


2353 - 


468 


0.25 - 


1.10 


Q 
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Table A- 3 — Continued 



IRDC core 


(K) 


T250 


/3 




L 


M 
{Mq) 






Classification^ 


G031.97+00.07 MM6 


10 - 39 


0.187 - 0.004 


2.5 - 


1.2 


52 - 1580 


2493 - 456 


1.58 - 


2.70 


Q 


G031.97+00.07 MM7 


8 - 46 


0.290 - 0.001 


2.5 - 


0.5 


10 - 347 


1142 - 91 


1.33 - 


3.37 


Q 


G033.69-00.01 MMIO 


10 - 41 


0.080 - 0.001 


2.5 - 


1.1 


26 - 827 


299 - 49 


0.64 - 


1.25 


Q 


G033.69-00.01 MMll 


11 - 39 


0.053 - 0.001 


2.5 - 


1.3 


33 - 987 


237 - 47 


0.73 - 


1.25 


Q 


G033.69-00.01 MM6 


16 - 43 


0.010 - 0.001 


1.6 - 


0.9 


58 - 690 


465 - 154 


< 0.01 - 


< 0.01 


Q 


G033.69-00.01 MM8 


11 - 39 


0.073 - 0.002 


2.5 - 


1.3 


35 - 1032 


1529 - 286 


0.09 - 


0.43 


Q 


G034.43+00.24 MM6 


7-48 


0.176 - 0.000 


2.4 - 


0.5 


2 - 93 


764 - 55 


0.04 - 


1.40 


Q 


G034.43+00.24 MM7 


16 - 43 


0.010 - 0.001 


1.5 - 


0.9 


15 - 187 


113 - 37 


0.74 - 


0.63 


Q 


G035.39-00.33 MM8 


11 - 49 


0.010 - 0.000 


1.1 - 


0.5 


2 - 49 


147 - 22 


< 0.01 


- 2.16 


Q 


G035.39-00.33 MM9 


18 - 38 


0.010 - 0.002 


1.6 - 


1.2 


17 - 173 


52 - 22 


< 0.01 - 


< 0.01 


Q 


G035.59-00.24 MM3 


13 - 48 


0.010 - 0.000 


1.4 - 


0.5 


4 - 52 


76 - 16 


< 0.01 - 


< 0.01 


I 


G048.65-00.29 MM2 


9 - 42 


0.133 - 0.001 


2.5 - 


0.9 


4 - 93 


131 - 18 


0.47- 


1.17 


Q 



denotes the classification based on tlie IRAC colors and the presence or absence of 24 /im emission. 



